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The  purpose  of  this  study  was  to  examine  the  various 
wavelength  dependent  issues  which  would  impact  the  selection 
of  an  operating  wavelength  for  a  ground  based  free  electron 
laser  weapon  system.  As  a  non-specialist  in  the  laser  field, 
my  approach  to  this  problem  was  to  develop  rather  simple 
scaling  relationships  and  combine  them  in  a  straightforward 
manner.  In  fact,  time  itself  precluded  a  more  thorough  treat¬ 
ment.  The  results  of  the  study  do  not  present  a  single  opti¬ 
mum  wavelength  but  show,  to  an  approximate  degree,  how  the 
selection  is  sensitive  to  various  parameters. 

No  analysis  is  a  solo  effort,  and  this  study  was  no  ex¬ 
ception.  My  thesis  advisor,  Lt  Col  John  Erkkila,  showed 
tremendous  patience  and  understanding  throughout  this  project 
and  deserves  my  highest  praise.  I  would  also  like  to  thank 
Dr.  Charles  Brau,  director  of  the  FEL  Project  Office  at  Los 
Alamos,  who  sponsored  this  research  and  provided  invaluable 
information  and  assistance.  Also,  much  thanks  is  due  Mr. 

Frank  Jinks  of  the  Wright  Aeronautical  Lab  who  helped  me  e- 
normously  in  getting  the  programs  FASCODE  and  LOWTRAN  up  and 
running.  Finally,  a  special  thank  you  goes  to  my  dearest 
friend,  India  Radford.  Her  love,  constant  support  and  en¬ 
couragement  have  truly  made  the  entire  effort  possible. 


Douglas  E.  Kohlhepp 
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A  computer  model  was  developed  to  examine  the  wavelength 
dependent  relationships  which  would  impact  in  a  selection  of 
the  operating  wavelength  for  a  ground  based  free  electron  la¬ 
ser  weapon  system  employed  in  a  strategic  role.  The  program 
considers  the  device  efficiency  of  the  free  electron  laser, 
the  propagation  of  laser  radiation  in  the  atmosphere,  and 
the  coupling  of  laser  energy  into  the  target  material.  The 
program  employs  modified  versions  of  the  atmospheric  trans¬ 
mission  codes  LOWTRAN  and  FASCODE. 

The  optimum  laser  wavelength  for  a  specific  selection  of 
laser  parameters  is  determined.  However,  equally  acceptable 
wavelengths  could  be  found  with  a  .2  micrometer  region  sur¬ 
rounding  this  optimum  value.  In  addition,  the  wavelength 
selection  process  is  very  sensitive  to  the  mean  sea  level 
altitude  of  the  laser  device  and  the  degree  to  which  atmo¬ 
spheric  turbulence  induced  beam  spread  can  be  compensated  for 
by  adaptive  optics. 


WAVELENGTH  SELECTION  FOR  A  GROUND  BASED 


FREE  ELECTRON  LASER  WEAPON  SYSTEM 


I.  Introduction 


Background 

Over  the  past  decade,  enthusiasm  over  the  military  uses 
of  directed  energy  weapons  has  increased  dramatically.  Spe¬ 
cifically,  considerable  interest  has  been  placed  on  the  de¬ 
velopment  of  ground  based  laser  weapon  systems.  Such  a  sys¬ 
tem  could  be  employed  as  an  anti-satellite  weapon,  or  as 
some  have  proposed,  in  a  strategic  defense  role  (1:17).  In 
this  latter  mission,  the  ground  based  laser  would  work  in  con¬ 
cert  with  high  quality  relay  mirrors  in  space  to  intercept 
ICBM's,  and  possibly,  high  altitude  bomber  aircraft. 

Basing  such  a  laser  weapon  on  the  ground  rather  than  in 
space  offers  many  advantages.  Generally,  size  and  weight  are 
not  the  crif '  ,al  considerations  on  the  ground  as  they  would  be 
in  space.  Ground  basing  allows  ease  of  servicing  and  repair. 
Perhaps  the  most  important  advantage,  however,  is  the  high 
survivability  that  basing  within  national  boundaries  offers 
over  space  basing. 

Yet,  the  ground  basing  of  a  laser  weapon  poses  a  very 
severe  disadvantage  as  well.  Such  a  weapon  system  would  be 
required  to  propagate  a  laser  beam  through  the  intervening 


atmosphere.  Atmospheric  attenuation,  turbulence,  and  self- 
induced  thermal  blooming  all  contribute  to  reduce  the  laser 
energy  that  reaches  the  target.  These  effects  are  very  much 
dependent  on  the  wavelength  of  the  laser  radiation;  slight 
changes  in  wavelength  can  drastically  alter  the  beam's  abil¬ 
ity  to  propagate  through  the  atmosphere.  In  addition,  the 
degree  to  which  a  laser's  energy  is  absorbed  by  a  target  al¬ 
so  determines  the  overall  effectiveness  of  a  laser  weapon 
system.  This  phenomenon  is  dependent  on  the  laser  wave¬ 
length  as  well.  Finally,  there  may  be  efficiencies  within 
the  laser  itself  which  are  functions  of  wavelength. 

Previously,  designers  of  laser  systems  for  defense  appli¬ 
cations  were  restricted  to  only  the  specific,  discrete  wave¬ 
lengths  produced  by  conventional  lasers.  With  the  advent  of 
the  free  electron  laser  (FEL),  however,  a  designer  would  now 
have  the  ability  to  select  the  wavelength  which  would  maximize 
the  effectiveness  of  the  weapon  system. 

. The  free  electron  laser  produces  coherent  electromagnet¬ 
ic  radiation  through  interactions  between  a  relativistic  e- 
lectron  beam,  a  spatially  periodic  magnetic  field,  and  an 
optical  radiation  field  propagating  in  the  same  direction  as 
the  electron  beam.  It  is  unusual  in  that,  unlike  a  conven¬ 
tional  laser,  it  does  not  relay  on  energy  transitions  between 
atomic  or  molecular  energy  levels.  Because  of  this  fact,  the 
operating  wavelength  of  the  free  electron  laser  is  continuous¬ 
ly  selectable.  In  theory,  any  wavelength  from  the  microwave 


region  to  the  x-ray  portion  of  the  spectrum  is  possible  (2). 
The  free  electron  laser  is  still  very  much  in  the  experimental 
stage  of  development,  but  the  eventual  scaling  to  the  high 
power  levels  required  of  a  laser  weapon  appears  promising 
(3:60). 

Review  of  the  current  literature  and  discussions  with  the 
Free  Electron  Laser  Project  Office  at  the  Los  Alamos  National 
Laboratory  indicate  that,  previously,  no  one  has  studied  at¬ 
mospheric  propagation,  target  absorption,  and  FEL  device  ef¬ 
ficiency  issues  in  an  overall  system  analysis  with  an  eye  to¬ 
wards  selecting  the  wavelength  which  best  accomplishes  a  spe¬ 
cific  mission.  This  analysis  is  necessary  if  one  is  ever  to 
take  advantage  of  this  unique  tunability  feature  of  the  FEL. 

Research  Objectives 

The  objective  of  this  study  is  to  examine  the  many  wave¬ 
length  dependent  relationships  which  affect  the  overall  ef¬ 
fectiveness  of  a  ground  based  free  electron  laser  weapon  sys¬ 
tem  and  to  select  that  wavelength  which  provides  the  maximum 
effectiveness  for  the  system.  Specifically,  this  objective 
will  include  the  following  subtasks. 

A.  To  determine  the  significant  wavelength  dependent 
issues  which  affect  weapon  system  efficiency. 

B.  To  develop  a  parametric  computer  model  which  eval¬ 


uates  the  overall  effectiveness  of  a  ground  based 
FEL  at  a  given  wavelength. 


C.  To  validate  the  model,  test  model  assumptions  and 


identify  inaccuracies  in  the  model 

D.  To  analyze  a  range  of  wavelengths  in  hopes  of  find¬ 
ing  a  single  wavelength  or  a  small  range  of  wave¬ 
lengths  which  maximizes  the  system  of  effectiveness 
for  the  input  parameters  chosen. 

E.  To  identify  the  key  issues  which  affect  the  wave¬ 
length  selection  process. 

Basic  Assumptions 

The  following  basic  assumptions  will  be  made  in  the  an¬ 
alysis  : 

A.  The  free  electron  laser  weapon  system  is  ground  based 
in  a  mid-latitude  region  of  the  United  States  and  in 

a  relatively  pollutant  free,  low  humidity  environment. 

B.  The  FEL  weapon  system  is  being  employed  in  a  strategic 
mission,  i.e.  anti-satellite,  ballistic  missile  de¬ 
fense,  anti-aircraft. 

C.  The  laser  can  be  used  in  a  direct  attack  mode  for  a 
target  in  space,  or  in  concert  with  a  100%  reflec¬ 
ting  spaced  based  relay  mirror.  It  is  assumed  that 
the  optimum  wavelength  would  be  nearly  the  same  in 
either  case. 

D.  The  free  electron  laser  produces  an  output  with  a 
long  pulse  length  which  can  be  assumed  to  be  con¬ 
tinuous  wave.  Hence,  the  device  peak  power  is  es- 


sentially  equal  to  the  average  power  (2). 


E.  The  target  is  composed  of  common  aerospace  materials 
such  as  aluminum  alloys. 

F.  The  kill  mechanism  of  the  weapon  system  is  a  thermal 
burnthrough  of  the  target's  outer  skin. 

G.  The  wavelength  dependence  in  the  reflectances  of  op¬ 
tical  mirrors  will  not  be  considered. 

H.  Efficiencies  associated  with  power  and  cooling  re¬ 
quirements  will  not  be  considered.  It  is  assumed 
these  are  insignificant  for  a  ground  based  system. 

I.  Pointing  and  tracking  capabilities  will  not  be  ad¬ 
dressed. 

J.  Further  assumptions  concerning  specific  areas  will 
be  discussed  during  the  chapters  pertaining  to  mod¬ 
el  development. 

Basic  Methodology 

The  brunt  of  the  analysis  effort  will  be  in  the  devel¬ 
opment  of  the  computer  model.  The  model  will  emphasize  a  sys¬ 
tem  wide  approach  to  the  problem  and  will  cover  the  following 
three  main  areas  of  an  interest  in  a  ground  based  free  elec¬ 
tron  laser  weapon  system:  FEL  device  efficiency,  atmospheric 
propagation  and  target  interactions. 

The  wavelength  dependence  in  the  device  efficiency  arises 
from  the  laser's  ability  to  trap  and  extract  energy  from  the 
relativistic  electron  beam  and  from  certain  restrictions  with- 


in  the  laser  cavity  (2).  Current  engineering  measures  which 
are  used  to  enhance  the  overall  FEL  device  efficiency  need  to 
be  considered  as  well. 

A  simplified  propagation  sub-model  will  be  developed  to 
address  the  propagation  issues.  It  will  consider  atmospher¬ 
ic  absorption  and  scattering  from  both  molecular  and  aerosol 

constituents  in  the  atmosphere.  In  addition,  the  spreading 
•  % 

of  the  laser  beam  due  to  diffraction,  jitter,  and  atmospher¬ 
ic  turbulence  will  be  included.  Finally,  the  model  will  ad¬ 
dress  the  reduction  of  the  laser  beam  irradiance  from  self- 
induced  thermal  blooming. 

Generally,  shorter  wavelengths  are  absorbed  by  metallic 
targets  more  efficiently  than  longer  wavelengths.  The  basic 
theory  behind  such  an  assumption  is  the  Drude  free  electron 
model,  and  thus,  will  form  the  foundation  for  this  portion  of 
the  analysis. 

Each  of  the  three  areas  just  discussed  can  be  described 
through  a  relationship  between  laser  wavelength  and  a  dimen¬ 
sionless  value  of  efficiency  between  zero  and  one.  The  pro¬ 
posed  computer  model  will  multiply  these  three  values  together 
to  form  an  overall  system  efficiency  as  described  below. 

n  system  *  n  FEL  x  n  propag  in,  x  11  target 

This  overall  efficiency  is  the  system  measure  of  effec¬ 
tiveness  we  wish  to  maximize  with  respect  to  laser  wavelength. 
Wavelength  selection  will  be  simply  based  on  the  maximum 


value  of  overall  system  efficiency. 


At  certain  times  in  the  analysis  process,  it  is  more  use¬ 
ful  to  employ  the  term  wavenumber  rather  than  wavelength. 

Both  are  spectral  representations  of  electromagnetic  radiation. 
This  report  will  use  the  terms  wavelength  and  wavenumber  al¬ 
most  interchangeably,  the  conversion  between  the  two  being 
governed  by  the  following  equation. 

4 

WAVENUMBER  (  on-1  )  =  A  /  ■  „  x  (2) 

\  /  WAVELENGTH  (micrometers)  v  ' 

Organization  of  the  Report 

In  the  following  three  chapters,  the  theory  supporting 
the  modeling  of  FEL  device,  propagation,  and  target  coupling 
efficiencies  is  presented  and  pertinent  wavelength  dependent 
relationships  are  developed.  Chapter  V  presents  the  program 
which  will  be  used  in  the  analysis.  The  actual  results  of 
the  analysis  will  be  presented  in  Chapter  VI.  Finally,  Chap¬ 
ter  VII  summarizes  with  the  conclusions  and  recommendations 
from  the  study. 


II.  Free  Electron  Laser  Efficiency  Model 


The  free  electron  laser  is  unique  among  laser  devices 
in  that  it  does  not  rely  on  the  transitions  between  discrete 
energy  levels  of  atoms  or  molecules  to  generate  stimulated 
emission.  Instead,  coherent  electromagnetic  radiation  is  pro¬ 
duced  through  the  interactions  of  a  relativistic  electron  beam, 
a  spatially  periodic  magnetic  field  or  'wiggler',  and  an  op¬ 
tical  radiation  field  propagating  in  the  direction  of  the  e- 
lectron  beam.  These  major  components  of  a  free  electron  la¬ 
ser  can  be  seen  in  a  cut-a-way  view  of  the  device  in  Figure 
1. 

Free  electron  lasers  have  only  been  in  existence  for  the 
past  seven  years.  Madey  and  others  at  Stanford  University 
first  demonstrated  applif ication  from  such  a  device  in  1976 
(4:717).  A  year  later,  the  same  group  demonstrated  laser  pow¬ 
er  from  a  cavity  resonator  operating  at  3.4  micrometers  (5:892). 
Since  those  early  breakthroughs,  much  has  been  published  in  an 
effort  to  advance  the  theory.  One  of  the  earlier  works  of  in¬ 
terest  was  that  of  Colson  (6)  who  described  free  electron  laser 
operation  totally  in  terms  of  classical  physics.  Also  of  note 
were  the  theoretical  proposals  to  enhance  the  efficiency  of 
the  device  through  novel  wiggler  designs  (7)  and  the  efforts 
to  optimize  certain  design  parameters  of  the  system  (8). 
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The  largely  qualitative  discussion  of  the  free  electron 
laser  which  follows  is  intended  to  provide  those  unfamiliar 
with  the  basic  concepts  of  the  device  an  understanding  of  the 
efficiency  relationships.  A  full  theoretical  treatment  is 
beyond  the  scope  of  this  study  and  the  interested  reader  is 
referred  to  the  many  excellent  references  on  the  subject. 

The  electron  beam  and  propagating  optical  beam  travel 
through  the  interior  of  the  wiggler  device.  The  wiggler  is 
characterized  by  a  magnetic  field  which  is  constant  in  magni¬ 
tude,  but  alternating  in  direction  down  the  length  of  the 

wiggler  with  some  spatial  period,  X  .  The  magnetic  field 

w 

vector,  B,  is  kept  perpendicular  to  the  path  of  the  electron 
beam  throughout.  This  is  generally  accomplished  through  a 
series  of  permanent  magnets  as  depicted  in  the  figure. 

As  the  electron  beam  travels  down  the  interior  of  the 
wiggler  and  encounters  the  magnetic  field  perpendicular  to 
its  path,  the  electrons  will  undergo  a  transverse  acceleration 
described  by  the  following  equation. 

3Y  (Ymv)  -  e  (v  x  I)  (3) 

where  m  «  electron  rest  mass,  Kg 

e  ■  electron  charge,  coulombs 
v  »  electron  velocity,  m/sec 
E  ■  the  magnetic  field  vector,  Tesla 
and  y  ■  the  relativistic  factor, 


Of  course,  the  direction  of  acceleration  due  to  the 
magnetic  force  is  normal  to  both  the  electrons'  path  and  the 
magnetic  field  because  of  the  cross  product  terms  in  the  equa¬ 
tion.  As  the  electrons  undergo  this  acceleration,  they  emit 
energy  in  the  form  of  electromagnetic  radiation.  The  result¬ 
ing  radiation  pattern  is  symmetric  about  some  mean  propaga¬ 
tion  vector  which  is  normal  to  the  direction  of  acceleration. 

In  this  case,  the  radiation  will  peak  in  the  direction  of 
the  electrons'  path.  Due  to  the  relativistic  nature  of  the 
electron  beam,  this  radiation  pattern  will  be  extremely  nar¬ 
row. 

If  the  magnetic  field  alternates  periodically  in  direction 
the  electrons  will  be  caused  to  accelerate  first  in  one  direc¬ 
tion,  then  another.  Therefore,  the  electrons  will  follow  an 
undulating  path  with  a  period  equal  to  the  spatial  period  of 
the  magnetic  field  (hence,  the  derivation  of  the  term,  wig- 
gler). 

Consider  the  radiation  produced  by  the  electrons  at  two 
points  along  their  path  as  depicted  in  Figure  2.  The  undulat¬ 
ing  motion  of  the  electrons  is  exaggerated  for  clarity. 
Actually,  the  electrons  deviate  very  little  from  their  orig¬ 
inal  path.  The  radiation  produced  at  Point  A  could  be  made 
to  combine  constructively  with  the  radiation  produced  at  Point 
B  with  some  judicious  selections  of  parameters.  The  time  nec¬ 
essary  for  a  photon  of  radiation  to  travel  the  distance  from 
A  to  B  is  simply 


where  c  is  the  speed  of  light.  An  electron,  traveling  close 
to  but  less  than  the  speed  of  light,  will  lag  behind  the 
photon.  Suppose  the  electrons  lag  a  distance  equal  to  ex¬ 
actly  one  wavelength  of  the  propagating  radiation.  In  this 
case,  the  radiation  they  produce  at  point  B  will  be  in  phase 
with  the  radiation  produced  at  point  A.  The  electrons'  time 
to  travel  from  A  to  B  is 


Te  -  Aw/v  (5) 

If  the  difference  between  the  two  times  corresponds  to  one 
wavelength  of  light,  A^  ,  we  have 


Aw  -  Aw  \  Al 
\  V  c  /  c 


(6) 


V 

c 


Aw 


(7) 


V  \ 

using  the  following  relativistic  relationship  and  its  Taylor 


series  expansion 


v 
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and  ignoring  the  insignificant  terms 
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Figure  2.  Electron/Wiggler  Field  Interactions 


we  arrive  through  simple  algebraic  operations  with  the  ex¬ 
pression 

X 

s  pf  do) 

Equation  (10)  defines  the  ’'resonance”  condition  for  a 
free  electron  laser  and  is  one  of  the  most  important  equa¬ 
tions  in  the  theory.  It  points  up  perhaps  the  most  remark¬ 
able  characteristic  of  the  free  electron  laser:  tunability. 
By  changing  either  the  wiggler  period,  A  ,  or  the  electron 
kinetic  energy,  y,  one  can  specify  the  wavelength  of  the 


laser  radiation.  Fox  example,  a  100  MEV  electron  beam  equates 
to  a  Y  of  about  196.  For  a  nominal  wiggler  period  of  5  cm,  a 
laser  output  in  the  visible  region  with  a  wavelength  of  .65 
micrometers  could  be  produced. 

This  simplistic  overview  is  intended  only  to  illustrate 
the  basic  operation  of  the  free  electron  laser.  In  reality, 
equation  (10)  is  applicable  only  for  weak  magnetic  fields 
where  one  can  assume  the  electron  velocity  is  not  perturbed 
appreciably  by  the  wiggler.  For  stronger  fields,  this  does 
not  hold,  and  the  resonance  condition  is  governed  by  the  fol¬ 
lowing  equation  (9:245). 
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where  Aw  is  the  dimensionless  wiggler  vector  potential  given 


by 


w 


2  me 
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Brau  (10:16)  employs  an  equivalent  expression  for  the 
wiggler  vector  potential  based  on  a  characteristic  length, 


^  ,  which  depends  on  the  permanent  magnet  material  and  wiggler 


construction,  and  the  separation  between  the  magnet  arrays, 
2h,  as  depicted  in  Figure  1. 


(13) 


Energy  Exchange  in  the  Wiggler.  Let  us  now  consider  the 


interaction  between  the  electrons  and  the  co-propagating 
electromagnetic  field  of  the  laser  beam.  This  interaction 
will  be  determined  by  the  Lorentz  force  equation. 

^  (Ymv)  =  e  (Er  +  v  x  Br)  (14) 

Here,  Er  and  Br  refer  to  the  electric  and  magnetic 
fields  of  the  laser  radiation.  The  magnetic  field  accomplish¬ 
es  no  work  on  the  electron;  therefore,  the  rate  of  change  of 
the  kinetic  energy  of  the  electron  can  be  described  by  the 
following  expression  (11:9). 

3t  (y)  =  3  ^  *  V  (15) 

me 

Equation  (15)  also  gives  rise  to  the  inter-related  con¬ 
cept  of  phase  between  an  electron  and  the  laser  field.  The 
sign  of  the  change  in  electron  engergy  depends  on  the  dot  pro¬ 
duct  of  the  electron  velocity  and  the  electric  field  of  the 
laser  beam,  and,  therefore,  on  the  spatial  positon  of  the 
electron  with  respect  to  the  light  wave.  Thus,  a  resonant 
phase,  in  this  simple  example  equal  to  zero,  can  be  defined 
(11:10). 

In  a  continuous  electron  beam  or  electron  bunch  of  large 
spatial  extent  compared  to  the  wavelength  of  the  light,  the 
phase  of  the  interaction  is  completely  random.  If  all  of  the 
electrons  are  at  the  resonance  energy,  approximately  half  will 
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Figure  3.  FEL  Power  and  Gain  vs.  Electron  Energy  (11:10) 


lose  energy  to  the  light  wave  and  half  will  gain  energy  from  the 


light  wave.  Hence,  at  resonance,  no  net  transfer  of  energy  takes 


place.  However,  if  the  electrons  are  injected  into  the  wiggler 


at  an  energy  slightly  higher  than  the  resonance  energy  for  a 


given  wavelength,  then  those  electrons  that  absorb  energy  from 


the  light  wave  are  pushed  to  higher  energies,  farther  from  the 


resonance  condition  and  proper  phase  with  the  light  wave.  Those 


electrons  that  lose  energy  to  the  wave  are  pulled  closer  to 


resonance  and  remain  in  phase  for  a  longer  period  of  time.  Thus, 


there  is  a  net  positive  energy  exchange  from  the  electron  beam 


to  the  laser  field. 


Figure  3  depicts  both  the  stimulated  emission  curve  and 


Figure  4.  FEL  Energy/Phase  Relationship  (7:99) 
the  small  signal  gain  curve  for  a  free  electron  laser.  As 
stated  above,  for  electron  beams  at  the  resonant  energy,  the 
laser  gain  is  zero.  For  higher  energies,  the  laser  gain  is 
positive  representing  transfer  of  energy  from  the  electron 
beam  to  the  laser  field. 

Those  electrons  that  are  decelerated  towards  the  resonant 
energy  are  eventually  captured  or  trapped  by  the  laser  field 
and  will  oscillate  about  the  resonant  values  of  energy  and 
phase.  Simply,  if  an  electron  is  decelerated  below  the  reson¬ 
ant  energy  and  phase,  it  will  now  begin  to  absorb  energy  from 
the  laser  field  as  described  by  equation  (15),  and  thus,  will 
be  accelerated  back  to  the  resonant  energy.  Therefore,  'buck¬ 
ets'  of  electrons,  trapped  by  the  laser  beam,  are  formed  with 
the  finite  dimensions  of  height  and  width.  Figure  4  depicts 
the  energy/phase  diagram  of  this  interaction.  The  area  with¬ 
in  the  closed  curves  represents  those  values  of  electron  en¬ 
ergy  and  phase  that  result  in  electrons  being  captured.  These 
closed  curves  effectively  define  the  height  and  width  of  the 


electron  buckets  as  well. 


Efficiency  enhancements 

Some  of  the  earliest  experiments  on  the  free  electron  la¬ 
ser  succeeded  in  extracting  less  than  .2%  of  the  electrons' 
energy  (4:720).  Since  those  times,  many  schemes  have  been 
proposed  to  enhance  the  efficiency  of  the  free  electron  laser. 
Perhaps  the  most  effective  has  been  the  development  of  the 
variable  parameter  wiggler. 

In  the  discussion  which  preceded,  it  was  assumed  the 
spatial  period  and  the  strength  of  the  magnetic  field  in  the 
wiggler  remained  constant.  Such  an  assumption  drastically 
limits  the  amount  of  energy  that  can  be  extracted  from  the 
electron  beam.  After  the  electrons  decelerate  in  the  wiggler 
and  become  trapped,  no  net  energy  transfer  between  the  elec¬ 
trons  and  laser  field  takes  place.  However,  if  the  value  of 
the  resonant  energy  is  slowly  (and  adiabatically)  decreased 
along  the  length  of  the  wiggler,  greater  amounts  of  energy  can 
be  extracted  from  the  beam  (7:108).  From  the  resonance  con¬ 
dition  defined  by  equation  (11),  it  can  be  seen  that,  holding 
the  laser  wavelength  constant,  the  resonant  energy  can  be  re¬ 
duced  by  decreasing  the  wiggler  period  and/or  decreasing  the 
wiggler  magnetic  field  potential. 

For  an  adiabatic  deceleration  of  the  electrons,  the  buck¬ 
et  size  and  shape  must  vary  slowly  along  the  length  of  the 
wiggler.  This  can  be  assured  if  the  resonant  phase  is  held 


constant.  Kroll  has  shown  that  this  resonant  phase,  ¥  , 
be  found  from  the  following  expression  (12:116). 
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where  z  is  the  position  along  the  wiggler  and  A  is  the  dimen 
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sionless  laser  field  equal  to 
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Here,  ZQ  is  the  impedance  of  free  space  =  (uq/£q)^  =  377  ft, 

r\ 

and  1^  is  the  laser  intensity  in  Watts/m.  Note  that  the  res 
onant  phase  reduces  to  a  value  of  zero  when  the  resonant  ener 
gy  is  constant  as  in  the  case  of  a  constant  parameter  wiggler 
In  the  case  of  variable  parameter  wigglers,  it  is  useful 
to  think  of  electrons  being  captured  in  buckets  and  then  the 
buckets  being  decelerated  to  extract  energy.  No  longer  is  it 
necessary  to  inject  the  electrons  at  an  energy  above  the  res¬ 
onant  energy.  Instead,  the  electrons  are  injected  at  the 
resonant  energy  and  are  immediately  captured.  These  buckets 
of  captured  electrons  are  then  slowly  decelerated  by  decreas¬ 
ing  the  resonant  energy  of  the  wiggler'.  With  such  a  model, 
one  can  define  a  capture  efficiency,  Hq,  and  a  deceleration 
efficiency,  of  the  free  electron  laser  device.  Then  the 
total  extraction  efficiency,  which  measures  the  device's 
ability  to  transfer  the  kinetic  energy  of  the  electrons  to 
the  radiation  field  of  the  laser,  is  simply  the  following 
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Note  that  the  extraction  efficiency  scales  inversely 
with  the  square  of  the  electron  energy,  indicating  the  greater 
difficulties  in  decelerating  higher  energy  beams.  Since,  the 
square  of  the  resonant  energy  is  likewise  inversely  propor¬ 
tional  to  the  laser  wavelength,  equation  (22)  shows  that 
higher  device  efficiencies  result  for  longer  wavelengths  as 
well . 


Cavity  Restrictions 

Of  course,  the  above  equation  represents  an  ideal  case 
where  there  exist  no  other  restrictions  to  the  design  of  the 
free  electron  laser  cavity.  Brau  and  others  at  the  FEL  Pro¬ 
ject  office  at  the  Los  Alamos  National  Laboratory  have  ex¬ 
tensively  studied  the  efficiency  of  the  free  electron  laser 
device  when  cavity  restrictions  such  as  diffraction  spreading 
of  the  laser  beam,  mirror  damage  and  mirror  alignment  are  con¬ 
sidered  (10:4-22).  The  rest  of  this  section  is  based  on  their 
work. 

Diffraction  Spreading.  The  maximum  length  of  the  wiggler 
magnet  and  thus  the  extraction  efficiency,  will  be  limited  by 
the  diffraction  spreading  of  the  laser  beam.  At  the  entrance 
and  exit  to  the  wiggler,  the  laser  field  must  not  impinge  on 
the  surfaces  of  the  premanent  magnets.  Brau  defines  the  fol- 
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lowing  quantities  (10:5). 
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(23) 


Here,  r  is  the  1/e  radius  of  the  beam  at  the  entrance 

(or  exit)  of  the  wiggler  and  r^  is  the  radius  at  the  beam 

waist.  It  is  assumed  the  beam  waist  falls  at  the  midpoint 

of  the  wiggler.  The  diffraction  spreading  is  governed  by  the 
•  \ 

following  equation 
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where  z  is  the  distance  from  the  beam  waist,  and  ZL  is  the 
Rayleigh  range, 


(25) 


Brau  (10:6)  defines  the  ratio  of  the  wiggler  aperture 
radius,  h,  and  the  beam  radius  at  the  entrance  to  the  wiggler, 
r,  as  simply 

n  -  £  (26) 


Finally,  the  electric  field,  equation  (17),  is  defined 
in  terms  of  the  laser  power,  P^,  and  a  characteristic  power, 
P0  (10:5). 


Here,  P0  equals  the  following. 
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Using  these  quantities,  Brau  determines  the  maximum  ex- 
ioi 

(10:6). 


traction  efficiency,  n  ,  based  on  a  maximum  wiggler  length 
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Mirror  Damage.  For  any  high  power  laser  device,  the 
maximum  power  that  can  be  generated  is  often  limited  by  the 
thermal  distortion  of  the  laser  cavity  mirrors.  This  is  also 
the  case  for  the  free  electron  laser.  For  an  assumed  Gaussian 
beam,  the  peak  laser  irradiance  at  the  mirror  surfaces  is 
given  by 


(30) 


where  r  is  the  laser  beam  radius  at  the  mirrors.  For  a  con- 
m 

tinuous  wave  laser  such  as  the  free  electron  laser,  mirror 
damage  will  most  likely  limit  the  time  averaged  irradiance  at 
the  mirror  surfaces  rather  than  the  peak  irradiance.  Since 
the  time  averaged  laser  power,  P^,  is  related  to  the  average 
electron  beam  current,  Tg,  by  (10:10) 
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Brau  shows  that  the  extraction  efficiency  of  the  device  is 
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further  restricted  by  mirror  damage  and  must  satisfy  the  fol¬ 
lowing  (10:11). 
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Here,  T  ,  is  the  limiting  average  laser  irradiance  on 
the  mirror  surfaces,  Ig,  is  the  peak  electron  beam  current, 
Lm  is  the  distance  separating  the  two  cavity  mirrors,  and  I 
is  a  characteristic  electron  beam  current  given  by 
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IQ  “  ffjj“  "  1356  (33) 

Cavity  Alignment.  It  can  be  seen  in  equation  (32),  that 
maximizing  the  extraction  efficiency  requires  g  >>  1.  or  the 
radius  of  the  laser  beam  at  the  wiggler  entrance  to  be  much 
greater  than  at  the  beam  waist.  This  requirement  necessates 
cavity  mirrors  with  short  focal  lengths,  thus  imposing  more 
severe  cavity  alignment  difficulties  (10:11). 

Figure  5  depicts  a  nearly  concentric  stable  laser  cavity 
with  mirror  radii  of  Rm«  A  small  angular  displacement,©  m, 
in  mirror  alignment  results  in  the  center  of  curvature  of  that 
mirror  being  displaced  by  an  amount  6  .  We  must  require  that 
the  tilt  of  the  laser  beam  axis  be  small  compared  with  the 
diffraction  angle  of  the  beam.  Making  the  approximation 
R  -  L  /2,  the  fractional  displacement  of  the  laser  spot  on 


The  capture  efficiency  of  the  wiggler  is  determined  by 
those  values  of  phase  at  which  electrons  are  captured.  Re¬ 


calling  Figure  4,  those  phase  values  will  fall  within  the  two 
extremes  that  define  the  bucket  width,  and  Therefore, 

assuming  every  phase  value  between  zero  and  2n  is  equally 
likely,  the  capture  efficiency  is  given  by  the  following  ex¬ 
pression  (10:1) 
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The  formula  for  the  decleration  efficiency  is  equally  simple 

(10:2). 
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Here,  Y  and Y r2  are  the  resonant  energies  at  the  entrance 
and  exit  of  the  wiggler.  Integrating  equation  (16)  with  re 
spect  to  z,  the  distance  along  the  wiggler,  yields 
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where  Lw  is  the  length  of  the  wiggler.  Therefore,  for  small 
amounts  of  deceleration,  the  extraction  efficiency  is  approx¬ 
imated  by  the  following  (10:4). 


the  mirrors  caused  by  mirror  misalignment  is  found  to  be 
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which  must  also  be  much  less  than  one  (10:11).  Imposing  this 
final  cavity  restriction  on  the  extraction  efficiency  of  the 
device,  leads  to  the  following  (10:11) 
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where  the  function  is  defined  by  the  expression 


FK(g,x) 
and  where 
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In  terms  of  the  variables  just  defined,  the  dimensionless 
wiggler  vector  potential  is  given  by  the  formula 
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and  the  resonant  electron  energy  at  the  entrance  of  the  wiggler 


by  the  expression 


1  +  A 


(40) 


Energy  Recovery 


Another  method  to  enhance  the  overall  efficiency  of  the 
•  \ 

free  electron  laser  is  to  recover  the  energy  of  the  electrons 
once  they  have  exited  the  wiggler.  This  energy  can  be  con¬ 
verted  to  radio  frequency  (rf)  power  which,  in  turn,  can  be 
used  to  accelerate  the  electrons  in  the  beam  prior  to  their 
entry  into  the  wiggler  (13:1).  In  general,  maximizing  the  ex¬ 
traction  efficiency  of  the  wiggler  does  not  maximize  the  over¬ 
all  efficiency  of  the  system  when  energy  recovery  techniques 
are  used  (13:1).  Therefore,  we  will  consider  this  last  fea¬ 
ture  of  the  device  in  order  to  develop  a  wavelength  dependent 
relationship  for  the  overall  system  efficiency. 

When  energy  recovery  is  used,  the  total  power  required  by 
the  electron  beam  accelerator  to  sustain  operation  of  the  de¬ 


vice  is  given  by 


P  A  +  PI  +  PM  +  nD  7B 


(41) 


where  PA  is  the  power  lost  in  the  accelerator  device,  P-j-  rep¬ 
resents  the  power  lost  in  the  injector,  ^  Pg,  is  the  power 
lost  in  the  wiggler  to  the  laser  field,  and  P^  is  the  power 
lost  during  energy  recovery  (13:1).  This  last  quantity  rep- 
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resents  the  fact  that  electrons  can  not  be  decelerated  below 


some  minimum  energy,  Yg,  before  control  of  the  beam  is  lost 
to  self-repulsive  effects.  Brau's  formulation  (13:2)  assumes 
the  power  lost  in  the  injector  is  proportional  to  the  power 
of  the  beam  in  the  wiggler;  therefore,  the  following  simplifi 
cation  is  used. 


5  +  p  =  YB 

M  1  Yrl 


B 


(42) 


Thus,  the  efficiency,  n  ,  for  converting  rf  power  to  laser 
power  is  given  by  the  following  formula  (13:3). 
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From  the  earlier  expression  for  the  extraction  efficiency, 
equation  (35),  the  overall  system  efficiency  becomes  (13:3) 
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where  the  variable  5  equals 
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Efficiency  Scaling 

In  the  design  of  a  free  electron  laser,  equation  (44) 
would  presumably  be  optimized  over  all  variables  subject  to 
the  laser  power,  P^,  being  fixed.  The  only  way  to  approach 


such  a  complex  optimization  problem  is  through  an  iterative  pro 
cedure  on  a  computer.  Fortunately  for  this  study,  this  proce¬ 
dure  has  already  been  accomplished. 

Brau  has  assumed  that  all  parameters  except  the  resonant 
phase,  ¥  ,  can  be  satisfactorily  optimized  by  first  maximizing 
the  extraction  efficiency  of  the  device  (13:5).  Having  done 
this,  equation  (44)  can  then  be  optimized  as  a  function  of  the 
resonant  phase,  ¥  .  The  process  proceeds  as  follows.  The 
parameter  <,  described  by  equation  (37),  is  used  as  input  to 
the  first  optimization,  from  which,  the  maximum  value  of  F ^  is 
determined  over  the  variables  g  and  x.  Figure  6  illustrates 
this  procedure,  where  for  a  fixed  <,  different  values  of  g  and 
x  are  iteratively  combined  to  fine  the  maximum  value  of  F<. 

Once  accomplished,  the  optimum  value  of  F<  is  used  to  compute 
5  .  With  6  so  fixed,  the  maximum  value  of  system  efficiency, 
nr,  is  found  through  numerical  iteration  of 

The  results  of  Brau's  optimizations  are  depicted  in  the 
next  four  figures.  The  optimum  values  of  g  and  x  are  found  in 
Figure  7  for  values  of  <  and  the  corresponding  maximum  values 
for  the  function  FK  are  depicted  in  Figure  8.  Figures  9  and 
10,  respectively,  show  the  optimum  value  for  the  resonant  phase 
1*  ,  and  a  scaled  value  of  optimum  system  efficiency  for  the 
fixed  parameter  <5.  Of  interest,  the  dashed  lines  in  these 
last  two  figures  represent  the  optimum  values  had  just  the  ex¬ 
traction  efficiency  of  the  device  been  maximized. 


Figure  7.  Optimum  Values  of  g  and  x  (10:15) 


There  exists  no  optimum  wavelength  in  the  free  electron 
laser;  however,  the  choice  of  laser  wavelength  will  indeed  af¬ 
fect  the  optimum  value  of  system  efficiency.  Figure  11  depicts 
this  relationship  between  wavenumber  and  the  optimum  system  ef¬ 
ficiency.  The  approach  taken  in  this  study  will  begin  by  assum¬ 
ing  that,  given  a  desired  wavelength  for  the  free  electron  laser 
output,  all  other  parameters  will  be  chosen  to  maximize  the  sys¬ 
tem  efficiency  with  energy  recovery,  n  .  Nominal  values  for 
certain  engineering  parameters  for  a  high  power  free  electron 
laser  device  (i.e.  permanent  magnet  sizes,  distance  between 
cavity  mirrors,  cavity  alignment  tolerances,  etc.)  were  provid¬ 
ed  by  the  FEL  project  office  at  Los  Alamos  (14:1-17).  With 
these  nominal  values,  an  optimum  system  efficiency  is  determined 
for  a  given  wavelength.  We  are  assuming  these  nominal  values 
would  not  change  significantly  over  the  range  of  wavelengths 
under  study. 

Instead  of  accomplishing  a  full  optimization  of  the  system 
efficiency  on  the  computer  for  each  wavelength,  optimum  values 
will  be  determined  through  the  curves  found  in  Figures  7  thru 
10.  The  data  from  these  curves  were  extracted  and  inserted  in¬ 
to  a  least-squares  polynomial  curve  fit  routine  to  determine 
the  relationship  between  variables.  The  optimum  values  of  g 
and  x  can  thus  be  found  for  a  given  <  by  the  following  expres¬ 
sions  . 

g  -  (-3.6495  k2)  +  (1.4410  <)  +  1.76241 


(46) 


x  =  (-8.8528  +  (4.0284  k)  +  .52725 


(47) 


from  which  the  optimum  values  for  F, ,  A^,  Yr- >  and  thus, 6  , 
are  determined.  Then  the  value  of  6  is  used  to  find  the  scaled 
system  efficiency  by 

(48) 

Scaled  Efficiency  =  (11.9642  52)  +  (-3.79196)  +  .55284 


The  scaled  efficiency  is  then  converted  to  an  actual  efficiency 
by  the  relationship 


n 

r 


_ 1  +  1A 

Scaled  Efficiency 


(49) 


This  curve  fit  approach  was  validated  against  typical  ef¬ 
ficiency  results  provided  by  the  FEL  Project  Office  which  were 
obtained  by  the  full  computer  optimization  (14:9).  It  was 
found  that,  near  the  center  of  the  curve  in  Figure  11,  the  re¬ 
sults  from  the  two  approaches  agreed  within  2%.  At  either  end 
of  the  curve,  however,  the  differences  grew  to  approximately 
4%.  It  is  felt  that  this  error  is  quite  acceptable  and  well 
within  the  assumed  errors  within  the  rest  of  the  model. 


III. 


Development  of  the  Propagation  Model 


The  vast  majority  of  the  overall  system  efficiency  in  a 
ground  based  laser  weapon  lies  in  the  ability  to  propagate 
the  laser  energy  through  the  atmosphere  to  the  target.  Fur- 
thur,  most  of  the  important  effects  are  strongly  wavelength 
dependent.  One  must  consider  linear  effects  such  as  the 
absorption  and  scattering  due  to  the  molecular  and  aerosol 
constituents  in  the  atmosphere,  and  the  spreading  of  the 
beam  due  to  atmospheric  turbulence,  diffraction,  and  plat¬ 
form  jitter.  For  most  power  levels  of  interest,  the  non¬ 
linear  thermal  blooming  effects  caused  by  the  self-induced 
heating  of  the  atmosphere  along  the  path  of  the  laser  beam 
must  also  be  considered. 

Many  computer  codes  have  been  developed  to  model  these 
propagation  effects  with  widely  varying  degrees  of  complex¬ 
ity.  Camphausen  developed  a  simplified  high  energy  laser 
propagation  code,  HELIM,  which  is  suitable  for  use  on  a  hand 
held  calculator  (15).  Lilly's  propagation  code  (16)  and  the 
program  GUTS  (Ground  Up  To  Space)  used  by  the  Air  Force  Weap 
ons  Laboratory  (17)  both  rely  on  Fourier  optics  for  their 
calculations.  Lying  somewhere  between  in  complexity 
is  the  code  developed  by  Peckham  and  Davis  (18).  Their  code 
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has  the  advantage  that  it  employs  relatively  simple  alge¬ 
braic  expressions  and  is  better  suited  for  the  system  ana¬ 
lyst.  The  propagation  model  we  will  develop  for  this  study 
will  be  similar  in  form  and  sophistication  to  Peckham  and 
Davis's  code. 

Overview 

For  this  study,  we  will  assume  a  laser  beam  profile  that 
approximates  an  infinite  Gaussian  distribution.  This  irradi- 
ance  profile  is  given  by 

I  -  Io  exp( -2r2/a02)  (50) 

where  I  *  the  peak  irradiance  at  beam  center 

O 

r  »  the  lateral  distance  from  beam  center 
and  a  »  the  l/ez  radius  of  the  beam,  within 

0  which  86.5%  of  the  beam  energy  is  found 

Therefore,  at  the  exit  aperture  of  the  laser,  the 
average  irradiance  across  the  beam  equals 

Power  ^  .865  Pfc  (51) 

*AV  *  Area  it  a0 

where  Pfc  is  the  total  average  power  contained  in  the 
beam.  After  propagating  a  distance  z  from  the  device,  the 


average  irradiance  is  given  by 


iAV(»)  *  -865  v'“Zlrel 


(52) 


where  a  *  the  attenuation  coefficient  of  the  atmosphere 
Irel  ■  the  power  reduction  factor  due  to  thermal 
blooming 

and  a  =  the  average  spot  size  of  the  laser  beam  at  a 
distance  z 

The  average  spot  size  can  be  approximated  by  combining 
in  a  root-mean-squared  technique  the  calculated  spot  sizes 
resluting  from  the  effects  of  diffraction,  turbulence  and 
platform  jitter  (19:1480)* 

a2  »  a2  +  a2  +  a2  (53) 

D  t  j 

To  define  a  propagation  efficiency  for  a  given  wavelength 
we  will  calculate  the  average  irradiance  incident  on  the  tar¬ 
get  and  divide  by  the  average  irradiance  produced  at  the  laser 
device . 


(54) 


The  above  constitutes  an  overview  of  the  basic  approach 
we  will  take  in  developing  the  propagation  model.  Each  of  the 


wavelength  dependent  factors  that  affect  laser  beam  propa¬ 
gation  will  be  discussed  in  greater  detail  in  the  following 
sections . 


Laser  Beam  Propagation  in  a  Vacuum 

Even  without  the  presence  of  an  intervening,  attenua¬ 
ting  medium,  a  laser  beam  will  spread  simply  by  diffraction 

\ 

due  to  the  wave  nature  of  light.  This  effect  can  be  reduced, 
with  restrictions,  through  optical  focusing  of  the  output 
laser  beam.  For  an  assumed  infinite  Gaussian  beam  profile, 
the  spot  size  of  the  beam  at  a  distance  z  from  the  laser  can 
be  found  through  wave  diffraction  theory.  Yariv  gives  the 
following  results  for  such  a  calculation  (20:33). 


a2(Z)  =  a2 

n  0 


(55) 


where  a  equals  the  radius  of  the  output  optics  of  the  laser 
0 

device  and  ^  is  the  propagating  wavelength.  For  an  unfocused 
beam,  the  minimum  beam  radius  or  beam  waist  is  found  at  the 
output  optics.  For  a  focused  beam,  however,  the  beam  waist 
occurs  at  a  point  between  the  laser  and  the  focal  length  of 
the  output  optics.  The  spot  size  at  the  waist  for  a  focused 


beam  is  also  given  by  Yariv  as 
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f  X/  TTa2) 

(56) 


where  f  equals  the  focal  length  of  the  output  optics  (20:38). 
Also,  the  focal  range,  2  ,  the  distance  from  the  laser  to  the 
beam  waist,  is  related  to  the  focal  length  by  the  following 
expression  (20:38). 
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(57) 


A  curious  result  of  the  above  equations  should  be  noted. 

If  the  focal  length  is  equal  to  the  Rayleigh  range,  Z  ,  de¬ 
li 

fined  as 


2 


L 


(58) 


the  spot  size  at  the  beam  waist  is  minimized  and  would  be 
equal  to 


In  addition,  the  focal  range  is  maximized.  In  this  case,  the 
beam  waist  is  located  exactly  halfway  between  the  laser  and 
the  focal  length.  If  one  attempts  to  focus  at  a  range  be¬ 
yond  the  Rayleigh  range,  the  spot  size  at  the  beam  waist 
increases  and  the  location  of  the  beam  waist  draws  nearer  the 
source.  In  this  case,  the  spot  size  at  the  target  is  greater 
than  if  no  focusing  were  employed. 

The  effects  of  focusing  and  diffraction  beam  spread  have 
been  combined  into  one  equation  by  Gebhardt  (19:1480),  Camp- 
hausen  (15:13)  and  others.  The  1/e2  radius  of  a  Gaussian 
laser  beam  at  a  distance  z  from  the  source  is  given  by 


The  term,  8  ,  in  equation  (60)  is  known  as  the  beam  qual¬ 
ity  factor.  Beam  quality  is  used  to  incorporate  imperfections 
in  the  phase  uniformity  of  the  beam  at  the  exit  of  the  laser 
device.  The  term  is  usually  defined  as  'times  diffraction 
limited'  at  some  empirically  derived  quantity  greater  than 


Laser  Beam  Attenuation 


Transmission  of  a  light  source  through  an  attenuating 
medium  is  defined  by  Beer's  Law  (21:19).  From  the  reference 
point  of  z=0  to  a  point  z  along  the  propagation  path,  the 
fraction  of  power  remaining  in  the  beam  is  given  by 

1(g)  (61) 

I(zl0)  =  exP<"a  z> 

The  attenuation  term,  “  ,  represents  a  combination  of 
many  different  loss  processes.  Included  in  this  term  are  the 
absorption  and  the  scattering  losses  from  both  molecular  and 
aerosol  constituents  in  the  atmosphere.  Each  loss  mechanism 
can  be  represented  by  its  own  attenuation  coefficient.  Addi¬ 
tion  of  each  term  yields  the  total  attenuation  coefficient 
for  the  atmosphere. 

(62) 

aTotal  =  ^molecular  absorption  +  01  aerosol  absorption 

+  “aerosol  scattering  +  “molecular  scattering 

where  for  each  individual  process 


a  (  A ,  Z)  =  n(2)  a  (  \) 


(63) 


Each  ot  is  dependent  upon,  n,  the  number  density  of 


absorbers  or  scatterers  the  beam  encounters  along  the  propa¬ 
gation  path  which,  of  course,  is  a  function  of  altitude  in 
the  atmosphere  (21:41).  The  wavelength  dependence  of  each 
process  is  contained  in  the  o  term  which  measures  the  cross- 
sectional  probability  of  an  interaction  occurring  (21:41). 
Each  of  the  major  attenuation  processes  represented  in  equa¬ 
tion  (62)  will  be  briefly  discussed. 

Molecular  absorption.  Molecular  absorption  is  the  di¬ 
rect  result  of  the  discrete  quantum  energy  states  of  the  mole 
cules  in  the  atmosphere.  When  a  photon  'collides'  with  an 
air  molecule,  and  the  energy  of  the  photon  equals  the  energy 
necessary  to  excite  the  molecule  from  a  lower  energy  state  to 
a  higher  energy  state,  the  photon  will  be  absorbed  by  the 
molecule.  Naturally,  since  radiative  energy  is  a  function  of 
photon  wavelength  (E=hA/c),  the  spectral  absorption  curve  is 
also  a  function  of  wavelength.  In  addition,  each  discrete 
absorption  wavelength  or  'line'  is  widened  by  collisional  and 
Doppler  broadening  processes  (21:43). 

At  wavelengths  longer  than  .75  micrometers,  atmospheric 
attenuation  is  most  strongly  influenced  by  the  contribution 
of  molecular  absorption.  The  structure  of  the  spectral  trans 
mission  curve  in  this  region  exhibits  enormous  complexity  and 
irregularity,  reflecting  the  discrete  energy  levels  of  the 
molecular  absorbers.  Figure  12  is  a  low  resolution  plot  of 


atmospheric  transmission  for  the  short  wavelength  infra-red 
and  visible  regions  of  the  spectrum  and  illustrates  this  com¬ 
plexity.  For  the  region  of  interest  in  this  study,  the  major 
atmospheric  absorbers  are  water,  carbon  dioxide,  oxygen,  and 
ozone  (22:146). 

Aerosols .  As  shorter  wavelengths  are  approached  in  the 
visible  spectrum,  the  contribution  of  aerosol  attenuation  be¬ 
gins  to  dominate  the  spectral  transmission  curve  (21:56). 

This  is  in  spite  of  the  fact  that  aerosol  number  densities 
are  generally  many  orders  of  magnitude  less  than  atmospheric 
molecules.  The  theory  of  aerosol  absorption  and  scattering, 
based  on  the  work  of  Gustav  Mie,  treats  atmospheric  aerosols 
as  small  dielectric  spheres  (21:57).  The  interaction  of  a 
light  ray  with  an  aerosol  can  cause  the  light  ray  to  be  re¬ 
flected  or  refracted  out  of  the  propagation  path  or  to  be 
internally  reflected  and  absorbed  by  the  aerosol.  The  absorb¬ 
ing  and  scattering  properties  of  aerosol  particles  is,  there¬ 
fore,  a  function  of  the  particle  radius,  the  wavelength  of 
the  propagation  radiation,  and  the  difference  in  the  refrac¬ 
tive  indices  of  the  aerosol  and  the  surrounding  air.  Although, 
Mie  theory  is  well  understood,  major  difficulties  arise  in 
application  due  to  the  extreme  altitude  and  geographical  vari¬ 
ations  in  the  sizes,  shapes,  and  dielectric  constants  of  the 
atmospheric  aerosols.  Generally,  however,  aerosol  attenuation 


is  a  slowly  varying  function  of  wavelength.  As  Figure  13 
illustrates,  shorter  wavelengths  are  somewhat  more  strongly 
attenuated  than  longer  wavelengths.  Mie  theory  is  applicable 
when  the  size  of  the  aerosol  is  comparable  to  the  wavelength 
of  the  propagating  radiation.  When  the  radius  of  the  aerosol 
becomes  much  larger  than  the  laser  wavelength  (as  in  heavy 
clouds,  fog,  and  rain)  attenuation  increases  dramatically  and 
becomes  wavelength  independent  (21:63).  Therefore,  this  type 
of  propagation  effect  will  not  be  considered  in  this  study. 
Conversely,  when  the  laser  wavelength  is  much  larger  than  the 
scattering  particle  (such  as  individual  molecules),  attenua¬ 
tion  becomes  a  strong  function  of  wavelength  (21:64).  This 
process  is  called  molecular  or  Rayleigh  scattering. 

Rayleigh  Scattering.  Rayleigh  scattering  results  from 
the  interaction  of  the  oscillating  electromagnetic  field  of 
the  laser  beam  and  the  individual  electrons  of  an  atmospheric 
molecule.  By  assuming  the  laser  wavelength  is  much  greater 
than  the  physical  dimensions  of  the  molecule,  spatial  varia¬ 
tions  of  the  electromagnetic  field  over  the  molecule  can  be 
ignored.  The  propagating  laser  field  induces  oscillations  in 
the  bound  molecular  electrons  which,  in  turn,  radiate  energy 
at  the  same  laser  frequency.  However,  not  all  of  this  energy 
is  radiated  in  the  direction  of  propagation,  and  thus,  power 
is  lost  from  the  laser  beam.  Using  this  approach  results  in 


v-vy> 


a  cross-sectional  probability  for  Rayleigh  scattering  which 
is  proportional  toA"^  (21:67).  This  strong  wavelength  depen¬ 
dence  can  be  seen  in  Figure  13.  At  wavelengths  above  1.0 
microns,  Rayleigh  scattering  is  only  a  minor  contributor  to 
laser  beam  attenuation;  however,  in  the  short  wavelength  por¬ 
tion  of  the  visible  spectrum,  it  can  be  seen  that  Rayleigh 
scattering  is  comparable  in  magnitude  to  total  aerosol  atten¬ 
uation  . 

Atmospheric  Transmission  Codes.  To  tackle  the  complex¬ 
ity  of  atmospheric  attenuation,  several  computer  codes  have 
been  developed.  Two  such  codes,  developed  by  the  Air  Force 
Geophysics  Laboratory,  which  were  available  for  this  research 
were  the  Fast  Atmospheric  Signature  Code  (FASCODE)  (24)  and 
the  Low  Resolution  Atmospheric  Transmission  Code  (LOWTRAN) 
(25). 

Both  codes  employ  the  concept  of  a  'layered'  atmosphere 
such  that  within  a  given  layer,  atmospheric  porperties  are 
assumed  constant.  The  altitude,  pressure,  temperature,  water 
vapor  density,  and  ozone  density  for  the  1962  U.S.  Standard 
atmosphere  and  five  seasonal  model  atmospheres  are  provided 
as  basic  input  data  to  both  LOWTRAN  and  FASCODE.  The  five 
supplemental  models  offer  the  user  a  selection  of  varying  at¬ 
mospheric  conditions.  They  include  models  which  represent: 

1)  Tropical  (15  N  Latitude),  2)  Midlatitude  Summer  (45  N, 
July),  3)  Midlatitude  Winter  (45  N,  January),  4)  Subartic 


Summer  (60  N,  July),  and  5)  Subartic  Winter  (60  N,  January). 
There  are,  however,  major  differences  between  these  two  codes. 

LOWTRAN  is  designed  to  provide  low  to  moderate  spectral 
resolution  results  to  the  user.  It  computes  atmospheric 
transmission  in  increments  of  five  wavenumbers  and  the  data 
for  each  wavenumber  has  been  averaged  over  twenty  wavenumbers 
(25:3).  It  does  consider  in  its  calculations,  all  loss  mech¬ 
anisms  which  comprise  the  attenuation  coefficient  of  equation 
(62).  FASCODE,  on  the  other  hand,  offers  the  highest  pos¬ 
sible  resolution  by  adding  the  contributions  from  individual 
molecular  absorption  lines  to  compute  transmission  at  a  given 
wavenumber;  however,  for  the  region  of  interest  in  this  study, 
it  does  not  include  aerosol  loss  mechanisms,  ozone  absorption, 
or  Rayleigh  scattering  (24:14)  Therefore,  neither  code  pro¬ 
vides  the  moderate  to  high  resolution  data  and  includes  every 
attenuation  process  which  is  needed  for  this  analysis.  In 
order  to  achieve  these  study  objectives,  it  will  be  necessary 
to  use  both  codes:  FASCODE  for  the  complex  molecular  absorp¬ 
tion  determination  and  LOWTRAN  for  the  more  slowly  varying 
functions  of  aerosol  absorption  and  scattering,  and  Rayleigh 
scattering.  The  data  generated  by  each  code  will  then  be  used 
in  the  model  developed  for  this  study.  Each  code  will  now  be 
discussed  in  greater  detail. 

LOWTRAN.  In  the  visible  region  of  the  spectrum,  aerosol 
absorption  and  scattering  become  very  important  contributors 


to  the  overall  attenuation  of  the  laser  beam.  LOWTRAN  employs 
an  aerosol  model  which  is  divided  into  four  regions,  each  con¬ 
taining  a  different  type  and  distribution  of  aerosols  (25:21). 
In  the  boundary  layer  region  from  sea  level  to  2  km,  the  size 
and  composition  of  the  aerosols  are  assumed  to  be  invariant 
with  altitude,  changing  only  as  a  function  of  relative  humid¬ 
ity  and  the  user  selected  environmental  conditions  (25:21). 

In  LOWTRAN,  the  user  may  chose  from  Rural,  Urban,  and  Mari¬ 
time  boundary  layer  aerosol  models,  and  further  specify  in 
kilometers  the  meteorological  range  or  visibility  at  the  sur¬ 
face.  The  number  density  of  aerosols  in  this  lowest  region 
is  assumed  to  vary  exponentially  with  altitude  but  only  under 
high  visibility  conditions.  Under  low  visibilities,  aerosol 
concentrations  become  independent  of  altitude  up  until  the 
altitude  of  the  first  temperature  inversion  (approximately 
1  km)  and  decrease  dramatically  above.  Above  this  boundary 
layer  in  the  upper  tropospheric  region  (2  to  10  km),  the  size 
and  distribution  of  aerosols  become  less  sensitive  to  surface 
conditions  and  seasonal  influences  dominate  (25:22).  In 
summer  months,  tropospheric  aerosol  concentrations  are  some¬ 
what  higher.  In  LOWTRAN' s  third  aerosol  region,  the  lower 
stratospheric  (10  to  30  km),  particle  size  and  distribution 
is  determined  by  the  degree  of  recent  volcanic  activity,  which 
at  times,  can  increase  aerosol  concentrations  by  a  factor  of 
100  over  the  normal  background  (25:34).  This  region  is  also 


scaled  by  LOWTRAN  to  compensate  for  seasonal  variations. 
Finally,  at  altitudes  above  30  km,  aerosols  composed  of  mete¬ 
oric  and  cometary  dust  form  the  major  component  of  the  normal 
upper  stratospheric  aerosol  model  in  LOWTRAN  (25:36).  Figure 
14  shows  the  vertical  profiles  of  the  aerosol  constituents 
used  in  LOWTRAN. 

The  cross-sectional  probabilities  for  aerosol  absorption 
and  scattering  are  computed  in  LOWTRAN  based  on  two  size  dis¬ 
tributions  of  aerosols,  the  so-called  accumulation  (small) 
and  coarse  (large)  particles  (25:25),  the  result  being  that 
the  absorption  and  scattering  (also  referred  to  as  extinction) 
coefficients  are  not  smooth  functions.  This  fact  can  be  seen 
in  Figure  15  which  depicts  the  coefficients  for  the  Rural 
aerosol  model,  normalized  to  1.0  at  .55  micrometers.  These 
coefficients  are  also  scaled  in  LOWTRAN  to  the  local  relative 
humidity.  As  the  relative  humidity  increases,  water  vapor 
condenses  out  of  the  atmosphere  onto  the  existing  aerosol 
particles.  This  condensed  water  thus  increases  both  the  size 
and  effective  refractive  index  of  the  aerosol  particles 
(25:26).  This  can  be  seen  as  well  in  Figure  15. 

LOWTRAN' s  treatment  of  molecular  scattering  is  fairly 
straightforward.  The  number  density  of  scatterers  is  scaled 
to  the  local  pressure  and  temperature  of  the  atmosphere.  The 
cross-sectional  probabilities  for  each  wavenumbers  were 
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Figure  14.  Aerosol  Concentrations  vs  Altitude 
in  LOWTRAN  (25:24) 


.  <*. 


obtained  from  least-square  curve  fits  of  molecular  scattering 
data  and  closely  follow  the  theoretical  functions.  The  errors 
in  this  curve  fit  are  less  than  .5%  in  the  region  of  interest 
of  this  study  (25:57). 

LOWTRAN  will  also  provide  molecular  absorption  data  for 
ozone  in  the  visible  and  ultraviolet  region  of  the  spectrum 
since  FASCODE  data  does  not  extend  into  this  region.  Figure 
16  depicts  the  two  spectral  regions  where  ozone  absorption  is 
found.  Of  note  is  the  fact  that  the  ozone  absorption  in  the 
ultraviolet  region  is  much  stronger  than  the  weak  ozone  ab¬ 
sorption  in  the  mid-visible.  Ozone  does  not  vary  exponential¬ 
ly  with  altitude  (25:48).  Instead,  ozone  concentrations  peak 
at  an  altitude  of  approximately  18  km,  the  so-called  ozone 
layer.  Each  of  the  six  atmospheric  models  in  LOWTRAN  contains 
separate  ozone  profiles.  For  the  ultraviolet  and  visible 
ozone  molecular  absorption  lines,  the  cross-sectional  proba¬ 
bilities  were  digitized  from  the  spectral  curves  of  McClatchey 
et  al,  in  intervals  of  500  wavenumbers  and  200  wavenumbers 
respectively  (25:48).  Therefore,  there  is  considerable 
interpolation  required  in  LOWTRAN  to  achieve  the  advertised 
five  wavenumber  resolution. 

Continuum  absorption  due  to  collision  induced  line 
broadening  of  water  vapor  and  nitrogen  is  available  in  both 
FASCODE  and  LOWTRAN;  however,  the  regions  where  continuum 
absorption  is  important  are  outside  the  region  of  interest 
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Figure  16.  Ozone  Attenuation  in  LOWTRAN  (25:53) 


in  this  analysis  (25:61) 

For  LOWTRAN,  each  of  the  six  atmospheric  models  is  divi 
ded  into  32  layers:  in  1  km  layers  from  sea  level  to  25  km, 
in  5  km  layers  from  25  km  to  50  km  with  the  final  two  layer 
boundary  at  70  km  and  100. 

LOWTRAN  provides  only  a  value  for  total  transmission 
along  a  selected  path  in  the  program  output.  In  order  to 
model  non-linear  effects  in  this  study,  it  was  necessary  to 
have  absorption  and  scattering  information  for  points  along 
the  beam's  path.  Therefore,  LOWTRAN  was  modified  to  extract 
this  data  at  each  of  the  32  altitude  layers  and  write  the  re 
suits  to  a  data  file  for  use  later.  The  specifics  of  this 
program  modification  are  included  in  Appendix  A  to  this  re- 
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FASCODE.  As  mentioned  earlier,  FASCODE  computes  molecu¬ 
lar  absorption  by  adding  contributions  from  individual  absorp¬ 
tion  lines.  At  lower  altitudes,  collision  broadening  mechan¬ 
isms  dominate  the  overall  width  of  the  molecular  absorption 
lineshape.  In  the  lower  atmospheric  pressures  of  higher  alti¬ 
tudes,  the  weaker  Doppler  broadening  mechanisms  will  begin  to 
dominate.  FASCODE  convolves  the  collision  broadened  line- 
shape,  approximated  by  a  Lorentz  function,  with  the  Doppler 
broadened  lineshape,  approximated  by  a  Gaussian  function,  into 
a  Voigt  lineshape  (24:10).  This  convolution  gives  an  effec¬ 
tive  description  of  the  overall  molecular  absorption  line- 
shape  at  any  altitude.  The  accuracy  of  this  approximation  is 
nominally  0.5%  with  a  maximum  error  of  3.0%  occurring  in  the 
wings  of  the  Voight  line  profile  (24:15).  This  maximum  error 
is  still  generally  less  than  the  errors  associated  with  the 
knowledge  of  both  the  lineshape  and  the  broadening  mechanism. 

For  program  efficiency,  FASCODE  samples  across  a  'panel' 
or  block  of  wavenumbers  of  interest  in  order  to  calculate 
atmospheric  transmission.  The  sampling  interval  is  determined 
by  the  average  width  of  the  molecular  absorption  lines  within 
the  block  of  wavenumbers  (24:11).  Since  the  average  linewidth 
changes  with  altitude,  the  sampling  interval  is  also  a  func¬ 
tion  of  altitude.  As  the  results  of  each  altitude  layer  are 
merged  with  preceding  layers,  FASCODE  employs  a  four-point 
Lagrangian  interpolation  to  match  the  different  sampling 
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intervals  (24:12) 


FASCODE  allows  the  user  to  either  specify  layer  bound¬ 
aries  or  employs  an  auto-layering  subroutine.  Of  course,  for 
this  study,  boundary  layers  were  chosen  to  coincide  with  the 
preset  layer  boundaries  in  LOWTRAN. 

Like  LOWTRAN,  FASCODE  had  to  be  modified  to  extract 
absorption  information  for  each  of  the  32  altitude  layers. 
These  modifications  are  described  in  Appendix  B.  In  one  re¬ 
spect,  this  was  a  simpler  process  since  FASCODE  computes 
absorption  data  one  layer  at  a  time.  However,  the  fact  that 
FASCODE' s  sampling  interval  varied  for  each  layer  required 
modifications  to  the  basic  program.  It  was  necessary  to  fix 
the  sampling  interval  for  each  layer  to  coincide  with  the 
desired  resolution  of  the  overall  analysis.  As  in  LOWTRAN, 
absorption  data  was  written  to  a  data  file.  In  addition,  the 
FASCODE  data  was  reconfigured  to  be  compatible  with  the  data 
format  of  LOWTRAN  and  the  main  program  used  in  this  study. 
Since  this  proved  to  be  a  non-trivial  matter,  this  additional 
program  is  presented  in  Appendix  C. 

Atmospheric  Turbulence 

Reduction  of  the  average  irradiance  on  the  target  will 
occur  due  to  turbulence  within  the  propagating  medium.  As  a 
laser  beam  propagates  through  the  atmosphere,  it  encounters 
numerous  cells  or  eddies  of  air  where  the  index  of  refraction 


is  slightly  different  from  that  of  the  neighboring  cells. 
These  small  changes  in  the  refractive  index  are  primarily  due 
to  the  small  variations  (.1  to  1.0  degrees  Kelvin)  in  the 
temperature  of  these  eddies  of  air  (21:70).  Irregular  solar 
heating  of  the  Earth's  surface  produces  these  temperature 
variations  in  the  atmosphere.  Naturally,  when  a  laser  beam 
encounters  a  cell  with  a  non-uniform  refractive  index,  its 
direction  and  phase  are  perturbed  by  a  small  amount.  This 
process  is  entirely  random,  but  by  dealing  with  statisti¬ 
cally  averaged  quantities,  one  can  successfully  approximate 
the  effects  of  turbulence.  While  the  refractive  index  varia¬ 
tions  are  quite  small,  long  propagation  distances  can  produce 
cumulative  effects  which  are  very  significant  (21:70).  A 
small  angular  deflection  near  the  laser  source  can  translate 
to  miss  distances  of  hundreds  of  meters  for  targets  in  geo¬ 
synchronous  orbits. 

Atmospheric  turbulence  effects  can  be  separated  into  two 
components:  instantaneous  small  scale  beam  spread  and  beam 

wander  (21:71).  When  the  lateral  diameter  of  the  refractive 
index  cells  are  small  compared  to  the  local  diameter  of  the 
beam,  turbulence  causes  small  scale  distortions  in  the  size 
and  shape  of  the  laser  beam  spotsize.  This  short  term  beam 
spread  component,  therefore,  describes  the  beam  spotsize  at 
any  instant  of  time.  Similarly,  when  the  size  of  the  turbu¬ 
lent  cell  is  large  compared  to  the  beam  diameter,  the  entire 


beam  follows  a  refracted  path.  This  component  of  turbulance 
thus  causes  the  beam  centroid  to  wander  about  some  central 
point  in  the  target  plane.  Each  component  is  typically  rep¬ 
resented  by  a  statistically  derived  characteristic  beam 
radius.  A  total  beam  radius  for  both  turbulence  effects  is 
then  found  by  a  root-sum-squaring  of  the  two  radii  (21:72). 

Refractive  Index  Structure  Parameter.  The  single  most 
important  factor  in  describing  turbulence  effects  in  the  at¬ 
mosphere  is  the  refractive  index  structure  parameter,  Cn2 
(26:20).  It  can  be  considered  a  measure  of  the  strength  of 
the  refractive  index  fluctuations  within  the  atmosphere. 

The  derivation  of  this  parameter  is  based  on  the  work  of 
Kolmogorov  and  Tatarskii  (26:20)  and  is  assumed  valid  over  a 
range  of  eddy  sizes--the  so-called  inertial  subrange  r  where: 

lo  <<  r  <<  L0  (64) 

Here  lo  and  L0  are  the  inner  and  outer  scales  of  turbu¬ 
lence  (note:  also  referred  to  as  the  microscale  and  macro¬ 
scale  of  turbulence  in  the  literature).  The  turbulent  eddies 
introduced  into  the  atmosphere  by  the  irregular  convective 
currents  from  the  Earth's  surface  and  wind  shear,  are  of  a 
large  scale.  Gradually,  the  turbulent  motion  of  these  cells 
causes  their  breakup  and  the  formation  of  smaller  and  smaller 
cells.  The  value  lo  corresponds  to  the  eddy  size  below  which 


dissipation  of  energy  in  the  eddy  through  viscous  effects 
becomes  important,  while  L0  is  the  largest  scale  size  for 
which  the  eddies  may  still  be  considered  isotropic.  Near 
ground  level,  10  may  be  on  the  order  of  millimeters  and  L0 
on  the  order  of  several  meters  (26:20). 

2 

As  one  might  expect,  the  value  of  Cn  decreases  with 
increasing  altitude  as  the  effects  of  surface  thermal  con¬ 
vection  and  frictional  wind  shear  become  less  noticeable. 
The  surface  of  the  Earth  can  have  such  a  strong  effect  that 
the  value  of  the  structure  parameter  near  the  surface  can 
easily  change  two  orders  of  magnitude  depending  on  the  so¬ 
lar  elevation  angle,  local  wind  conditions  and  cloud  cover 
(26:21).  Therefore,  predictive  models  of  the  behavior  of 
the  structure  parameter  with  altitude  at  best  represent 
nominal  values.  In  addition,  models  which  are  useful  near 
the  surface  are  of  little  applicability  in  the  so-called 
free  atmosphere. 

Near  the  surface,  the  structure  parameter  decreases 
proportionally  with  altitude  to  the  -2/3  power  during  quiet 
periods  in  the  atmosphere,  i.e.  dawn  and  dusk.  However, 
with  the  unstable  conditions  of  a  sunny  summer  day  at  noon¬ 
time,  a  more  appropriate  scaling  would  be  with  altitude  to 
the  -4/3  power  (26:22-25).  Wyngaard  et  al  have  developed 
such  a  model  for  the  latter  case  (27:1646-1650).  The 
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model  is  valid  from  the  surface  to  the  first  temperature 
inversion  in  the  atmosphere,  which  on  a  warm  summer  day  may 


be  over  1000  meters  high. 


(65) 


Cn2(h)  =  (dn/dT)2  ( 2xl0~3 )  Q4/3h~4/3 

14000  U.,.3  ~|  2/3 
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where  h 
dn/dT 


and  U , 


height  in  meters  above  the  ground 
the  change  in  refractive  index  per 
unit  change  in  temperature,  K_l  2 

the  upward  convective  heat  flux,  W/rn 
a  characteristic  frictional  velocity 
with  the  surface,  m/sec 


Typically,  U^.  is  an  order  of  magnitude  smaller  than 
the  local  wind  speed  and  may  be  approximated  by 


J*  -  0.35  ^h 


(dU/dh) 


(66) 


where  dU/dh  equals  the  vertical  gradient  of  the  mean  hori¬ 
zontal  wind  (28:6.11).  For  many  conditions,  the  term 
£h  (dU/dh)j  becomes  independent  of  altitude.  A  nominal 
value  for  this  factor  from  Wyngaard's  experimental  data 
over  the  plains  of  Kansas  would  be  1.0  m/sec  (27:1648). 

The  value  of  Q  depends  on  the  amount  of  sunlight  falling  on 
the  ground  and  the  properties  of  the  ground  itself.  As  an 
estimate  of  Q,  one  may  take 


Q  ■  Qo  Sin  5  -50 


(67) 
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where  €  is  the  solar  elevation  angle  and  Q0  is  a  constant 
containing  the  effects  of  cloud  cover  and  surface  composi¬ 
tion  (28:6.11).  For  a  clear  sky  and  dry  fields,  Qo  has  a 
value  of  400  w/m^  (28:6.12). 

An  approximate  formula  for  the  refractive  index  of  dry 
air  at  optical  frequencies  was  obtained  from  Clifford 
(26:10). 


1  +  (7.76xlO~5)P  +  (5.84xlO"I9)P 


(68) 


where  P  is  atmospheric  pressure  in  millibars,  T  is  temper¬ 
ature  in  degrees  Kelvin  and  *  is  the  wavelength  of  the  pro¬ 
pagating  light  in  meters.  Differentiating  equation  (68) 
with  respect  to  temperature  yeilds  an  appropriate  expres¬ 
sion  for  the  change  in  refractive  index  for  a  unit  change 
in  temperature. 


dn  _  (7.76xlO~5)P  ( 5 . 84x10~I9)P 


(69) 


At  altitudes  near  the  first  temperature  inversion, 
there  is  typically  a  local  peak  in  the  structure  parameter 
curve  (28:6.12).  However,  even  at  much  higher  altitudes 
in  the  free  atmosphere,  the  behavior  of  the  structure 
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parameter  is  characterized  by  sharp  peaks  of  possibly  two 
orders  of  magnitude  (28:6.18).  A  particularly  significant 
peak  occurs  in  the  9-15  kilometer  region  near  the  tropo- 
pause.  Yet,  beyond  the  tropopause,  the  structure  parameter 
quickly  approaches  zero  (28:6.19). 

Hufnagel  (26:24)  has  incorporated  this  highly  active 
tropopause  region  in  the  derivation  if  his  latest  model  of 
the  behavior  of  the  structure  parameter  in  the  free  atms- 
sphere.  The  model  is  valid  for  altitudes  above  the  first 
temperature  inversion  to  25  km.  The  mean  values  of  Cn  as 
a  function  of  altitude,  h,  is  given  by 

•  (70) 

Cn2(h)  *  2.7  <  (2.2xl0'53h10)  exp(-h/1000)  1 

+  10-16exp( -h/1500)> 

A  plot  of  Hufnagel' s  model  is  shown  in  Figure  17. 

As  stated  before,  the  region  which  most  influences 
turbulence  effects  in  the  far  field  is  near  the  surface 
where  the  laser  source  is  located.  Therefore,  this  study 
will  use  the  two  models  just  described  to  determine  the  be¬ 
havior  of  the  structure  parameter  with  altitude.  Wyngaard's 
model  will  be  used  from  the  surface  to  1000  meters,  an 
appropriate  value  for  the  first  temperature  inversion  for 
sunny  daytime  conditions,  above  which,  Hufnagel 's  model  is 
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employed  to  an  altitude  of  25  km.  Above  this  altitude, 

the  structure  parameter  can  be  assumed  to  be  zero.  The  two 

- 16  -2/3 

models  merge  at  a  value  of  approximately  1.5x10"  m  '  at 
an  altitude  of  1000  meters  when  a  mid-afternoon  sun  angle 
of  40  degrees  is  used. 

The  use  of  a  structure  parameter  model  near  the  sur¬ 
face  which  essentially  describes  a  worse  case  scenario  in 

\ 

terms  of  turbulence  effects  is  very  conservative  but  justi¬ 
fied  by  the  fact  that  the  operational  use  of  a  ground  based 
laser  weapon  system  would  require  an  immediate  response 
even  under  such  unstable  atmospheric  conditions. 

Turbulence  Induced  Beam  Radius.  The  theoretical 
development  of  the  effects  of  turbulence  on  a  propagating 
laser  beam  follows  the  work  of  Yura  (29).  The  key  parame¬ 
ter  in  his  formulation  is  P  ,  the  lateral  coherence  length. 

0 

This  quantity  can  be  thought  of  as  the  separation  between 
points  in  a  plane  perpendicular  to  the  propagation  direc¬ 
tion  where  a  spherical  wavefront  remains  coherent.  Yura 
presents  the  following  formula  for  this  quantity  valid  over 
the  propagation  range  Z^  <<Z<<Z2  (29:2771). 


,45  k  sec 
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(71) 
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where 


Z1  =  (.4  k2  Cn2  (Lo/2tt  )5/2)'1 


(72) 


Z2  -  (.4  k2  Cn^o5^3)-1 


(73) 


Here,  k  is  the  propagating  wavenumber  =  2  it/ A,  <p 

2 

is  the  zenith  angle,  and  Cn  is  the  refractive  index 

-2/3 

structure  parameter  in  units  of  meters  .  Note  that 

the  integration  of  the  structure  parameter  is  weighted 

more  heavily  toward  values  located  near  the  source  of  the 

laser.  Consider  visible  frequencies  and  a  constant  ground 

-  1 

level  value  for  the  structure  parameter  of  2.0x10  .  The 

propagation  ranges  of  validity  according  to  equations  (72) 

and  (73)  lie  between  .5  and  400,000  meters.  Propagating 

vertically  through  the  atmosphere  would  mean  average  value 
2 

of  Cn  several  orders  of  magnitude  below  a  ground  level 
value.  Therefore,  it  can  be  seen  that  Yura's  theory  is 
valid  over  the  ranges  of  interest  in  this  study,  and  in 
fact,  over  most  ranges  for  any  high  powered  laser  system 
study . 

Yura  assumes  an  infinite  Gaussian  profile  for  the 
laser  source  in  his  formulation.  Through  his  statistical 
treatment  of  turbulence,  the  beam  profile  in  the  far  field 


also  approximates  a  Gaussian  profile.  Yura  proceeds  by 
approximating  a  1/e  characteristic  total  beam  radius  de¬ 
pendent  on  the  lateral  coherence  length  parameter  (29:2771) 
Converting  this  formula  to  a  1/e  radius  yields  a  total 
beam  spread  radius  due  to  turbulence  effects  which  will 
be  used  in  the  model  for  this  study. 


The  effects  of  turbulence  might  be  compensated  for 
through  the  use  of  adaptive  optical  systems.  Without  cor¬ 
rection,  the  beam  characteristics  of  a  ground  based  laser 
device  would  quickly  be  dominated  by  turbulence  effects 
and  significantly  reduce  the  weapon's  utility.  The  degree 
to  which  this  compensation  can  be  accomplished  remains  an 
open  question  and  will  be  treated  parametrically  in  the 
model  developed  for  this  study. 

Effects  of  Jitter 

Vibrations  in  or  near  a  high  energy  laser  device  can 
cause  small  angular  deflections  in  the  laser  beam  at  the 
source.  Over  long  propagation  distances,  these  small  de¬ 
flections  can  translate  to  significant  errors  in  the 
target  plane.  These  vibrations  can  be  the  result  of 


frictional  forces  within  moving  parts  of  mechanical  de¬ 
vices  or  merely  the  relatively  constant  seismic  activity 
within  the  Earth.  Most  can  be  isolated;  however,  even  a 
free  electron  laser  device,  with  few  mechanical  parts, 
will  experience  a  minor  amount  of  jitter. 

These  jitter  effects  on  the  laser  beam  are  normally 
treated  statistically  in  laser  system  studies.  A  mean 
squared  radial  displacement  of  the  laser  spotsize  at  a 
propagation  distance  z  can  be  determined  by  the  following 
formula  (19:1480). 


a2  -  2  <V\  Z2 

J  '  x/ 


(75) 


where 

angle. 


(e2\ 

\  X/ 

This 


is  the  variance  of  the  single  axis  jitter 
formula  assumes  an  isotropic  jitter,  i.e. 


^9^  =  ^9^  .  Even  though  jitter  is  not  a  wavelength 


x'  '  y 

dependent  phenomenon,  it  will  be  included  in  this  study  in 
order  to  evaluate  non-linear  effects  through  the  atmosphere 


Thermal  Blooming 

When  a  high  power  laser  beam  passes  through  an  ab¬ 


sorbing  medium,  the  phenomenon  of  thermal  blooming  will  act 
to  distort  the  laser  beam  profile  and  could  significantly 


reduce  the  laser  beam  power  which  reaches  the  target.  In 
fact,  this  process  is  often  the  limiting  factor  in  deter¬ 
mining  the  maximum  amount  of  power  which  can  be  transmitted 
through  the  atmosphere  at  a  given  wavelength. 

The  essential  features  of  thermal  blooming  can  be 
explained  through  relatively  simple  physical  processes. 
Consider  an  ideal  gas  which  is  initially  in  thermal  equi¬ 
librium  and  with  an  ambient  index  of  refraction.  If  the 
gas  absorbs  radiation  at  the  wavelength  of  the  laser,  this 
absorbed  energy  gives  rise  to  local  heating  of  the  gas. 

The  increase  in  temperature  produces  a  small  increase  in 
pressure.  The  medium,  expanding  at  the  speed  of  sound  in 
response  to  this  pressure  imbalance,  decreases  in  density, 
and  thus,  causes  a  proportional  decrease  in  the  local  re¬ 
fractive  index. 

Consider  also  a  laser  beam  which  initially  has  a 
Gaussian  intensity  profile  as  in  Figure  18a.  On  the  axis 
of  the  beam,  heating  of  the  medium  will  be  at  a  maximum; 
therefore,  the  local  index  of  refraction  will  be  at  a 
minimum  value  as  shown  in  Figure  18b.  The  index  of  refrac¬ 
tion  returns  to  the  ambient  value  near  the  edges  of  the 
laser  beam.  Thus,  light  rays  near  the  center  of  the  beam 
will  be  refracted  radially  outward  toward  the  more  dense 
regions  of  the  gas.  The  laser  beam,  therefore,  diverges 
or  'blooms'  as  the  heating  of  the  gas  continues. 


An  interesting  complication  occurs  when  the  absorbing 
medium  is  flowing  transverse  to  the  direction  of  propaga¬ 
tion  by  means  of  either  a  natural  crosswind  or  a  slewing  of 
the  laser  device.  Under  these  conditions,  the  refractive 
index  profile  will  be  skewed  as  in  Figure  18c.  The  trans¬ 
verse  motion  of  the  medium  continuously  replaces  that  por¬ 
tion  of  the  medium  which  has  been  heated  by  the  laser  beam 
with  unheated  medium.  Therefore,  the  upwind  region  of  the 
beam  is  cooler  and  more  dense  than  the  downwind  region 
which  has  undergone  laser  heating  longer.  As  before,  light 
rays  are  refracted  toward  these  cooler  regions.  In  this 
case,  the  process  produces  an  asymmetric,  crescent  shaped 
beam  profile  depicted  in  Figure  18d. 

If  the  laser  beam  power  and  absorbtivity  of  the  atmo¬ 
sphere  are  sufficiently  large,  the  beam  will  simply  be 
obliterated  by  thermal  blooming  effects  after  propagating 
only  a  short  distance.  Studies  have  shown  that  this  can 
occur  on  time  frames  of  milliseconds  (30:222).  Yet,  for 
modest  power  levels  and  absorbtivity,  convective  cooling 
in  the  beam  and  the  transverse  flow  of  air  across  the  beam 
can  combine  to  create  a  steady-state  condition  for  thermal 
blooming  (31:63). 

Not  only  is  the  beam  expanded  by  the  self-induced 
thermal  blooming,  but  the  effective  center  of  the  beam  is 
displaced  in  the  direction  of  the  wind  (31:64).  This 
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effect  could  complicate  aiming  of  the  beam;  however,  for 
laser  beams  which  survive  the  effects  of  thermal  blooming, 
thermal  'bending'  effects  are  small  and  will  not  be  con¬ 
sidered  in  this  analysis. 

A  quantitative  description  of  the  thermal  blooming 
process  can  follow  several  approaches.  The  problem  can  be 
treated  through  sophisticated  wave  optics  -alculation ; 
however,  this  methodology  has  one  key  disadvantage.  Most 
non-specialists  would  quickly  lose  sight  of  the  approach 
and  the  basic  physical  processes  involved.  An  alternate 
avenue  might  be  the  one  taken  by  Gebhardt  and  Smith  (31). 
They  employed  relatively  simple  geometric  optics  approxi¬ 
mations  to  develop  a  thermal  blooming  'distortion'  parame¬ 
ter.  This  parameter  is  then  used  as  an  argument  to  enter 
a  series  of  curves  to  obtain  the  amount  of  irradiance 
reduction  on  the  target  due  to  thermal  blooming.  The 
curves  were  fit  from  experimental  data  and  results  from  de¬ 
tailed  wave  optics  computer  codes.  The  approach  to  thermal 
blooming  by  Gebhardt  and  Smith  is  the  one  which  will  be  de¬ 
veloped  here  and  used  in  the  model  for  this  study. 

As  laser  energy  is  absorbed  by  the  medium,  it  is  trans¬ 
mitted  to  cooler  surroundings  through  either  convection  or 
conduction.  The  change  in  temperature  across  the  medium  is 
governed  by  the  following  energy  balance  equation  (31:64). 
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where  P  =  atmospheric  density  in  kg/m 
0^  =  specific  heat  in  joules/kg-K 

V  =  wind  velocity  transverse  to  the 
propagation  path  in  m/sec 
a  =  absorption  coefficient  of  the  atmosphere 
K  =  the  conductivity  of  the  atmosphere 
and  I  =  the  laser  irradiance  in  watts/m2 

In  this  theoretical  development,  convective  heat 
transfer  is  assumed  to  dominate  over  thermal  conduction; 
therefore,  the  second  term  in  equation  (76)  is  ignored.  F 
thermore,  it  is  assumed  that  the  time  for  the  conversion 
of  absorbed  laser  energy  to  thermal  motion  is  small  com¬ 
pared  to  some  characteristic  convective  heat  transfer  time 
equal  to  a/Vx  where  a  is  the  local  radius  of  the  beam. 
Finally,  energy  losses  due  to  re-radiation  of  the  absorb¬ 
ed  laser  energy  are  neglected. 

The  refractive  index  of  the  medium  will  vary  about 


some  ambient  value,  n0,  due  to  the  small  spatially  vary¬ 
ing  temperature  changes  across  the  beam  (31:64).  For  an 
arbitrary  position  vector,  r,  we  have 


If  we  assume  these  variations  in  the  refractive  index 
are  small  and  occur  over  dimensions  large  compared  to  the 
laser  wavelength  (easily  satisfied  at  optical  frequencies), 
the  following  scalar  wave  equation  can  be  used  to  describe 
laser  propagation  (31:64). 


( V  2  +  k2n2  +  ik  a  )  U(r)  =  0 


(78) 


where  k  is  the  vacuum  wavenumber  and  U(r)  is  the  complex 
electric  field  of  the  laser  radiation. 

Geometric  optics  assumes  that,  for  sufficiently  short 
wavelengths,  one  can  describe  surfaces  of  constant  phases 
called  phase  fronts  such  that  normals  to  these  surfaces 
represent  the  direction  of  local  energy  flow  (32:235). 
Consider  as  a  trial  form  of  the  fixed  time  complex  electric 
field  of  the  laser, 


U(r) 


A(r)  exp 


(79) 


where  r  is  the  position  vector  of  a  point  along  a  given 
ray  of  light  and  S(r),  independent  of  the  magnitude,  A(r), 
is  a  scalar  function  of  position  which  locates  these  sur¬ 
faces  of  constant  phase  (32:235).  Substituting  equation 
(79)  into  Maxwell's  equations  yields,  in  the  short  wave¬ 
length  limit, 
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where  n(7)  becomes  the  position  dependent  value  of  refrac¬ 
tive  index  along  the  ray.  The  function  S(r)  is  often  call¬ 
ed  the  eikonal  and  equation  (80),  the  eikonal  equation. 

Both  n(r)  and  the  gradient  of  S(r)  are  normal  to  the  sur¬ 
face  of  the  phase  front.  From  Born  &  Wolf  (33:121-122), 
the  eikonal  equation  can  be  developed  into  a  differential 
equation  for  the  ray  vector. 


(81) 
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(  yS(r))  =  v  n(r) 


Here,  s  represents  the  scalar  measure  of  distance  along 
the  propagating  ray.  Note  that  a  constant  refractive  in¬ 
dex  implies  an  unvarying  gradient  of  S(r)  and  no  beam 
divergence . 

By  assuming  a  plane  wave  form  for  S(r)  and  the  parax¬ 
ial  ray  approximation,  i.e.  the  ray  vector  is  nearly  paral¬ 
lel  to  the  axis  of  propagation,  Gebhardt  and  Smith  developed 
an  irradiance  distribution  of  the  thermally  bloomed  laser 
beam  (31:65). 


TO'W,- '."A.1  T \  K TT '  V-'  -  ■  •  ■  V ■  '  ■  ■ 


I  (x,y,Z)  = 

Iu(x»y>z)e"  “  Z  exp 


=  Iu(x,y,Z)e-aZ  e" 

Where  I  is  the  unbloomed  irradiance  distribution  at  pro¬ 
pagation  distance  z  and  V  is  the  transverse  gradient 
operator  =  d/dx  x  +  d/dy  y.  The  attenuation  factor,  e~a* 
accounts  for  absorption  in  the  medium.  This  expression  con¬ 
veniently  describes  the  effect  of  refractive  index  varia¬ 
tions  across  the  laser  beam  in  terms  of  an  exponential  fac- 
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tor,  e  ,  that  modifies  the  unbloomed  irradiance,  Iu>  For 
a  Gaussian  beam  profile  given  previously  by  equation  (50), 
the  total  useable  power  within  the  1/e  radius  of  the  beam 
is  found  through  integration. 

(83) 
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therefore  one  can  write 
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The  change  in  temperature  across  the  beam  is  obtained  from 
equation  (76). 
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Combining  equations  (77)  and  (85)  yields  the  following 
(31:66). 
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where 
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Hence,  all  the  elements  in  equation  (82)  have  been 
defined.  Through  rather  straightforward,  albeit  laborious, 
mathematics  and  again  invoking  the  paraxial  ray  assumption 
and  ignoring  the  z  dependence  on  the  transverse  coordinates, 
x  and  y,  an  expression  for  V  can  be  determined  (18:32). 

(88) 

'F  =  NI  (^constants  independent  of  Z  ) 


Here,  NI  is  Gebhardt  and  Smith's  thermal  blooming  dis¬ 
tortion  parameter  defined  as 
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They  argue  that  the  reduction  of  laser  beam  irradiance 
due  to  thermal  blooming  effects  can  be  determined  solely  by 
this  single  distortion  parameter,  NI  and  have  developed 
an  empirical  model  based  on  this  parameter.  Figure  19 
shows  the  dependence  of  an  irradiance  reduction  factor, 

Irel,  on  the  distortion  parameter,  NI.  The  curve  was  de¬ 
veloped  from  experimental  data  points,  depicted  by  open 
circles  in  the  figure,  and  through  theoretical  results  of 
nonlinear  wave  optics  propagation  codes,  shown  as  solid 
circles.  The  curve  fit  results  predict  the  following 
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Figure  19.  Beam  Reduction  Factor,  Irel, 

vs.  Distortion  Parameter  (19:1484) 

The  elements  under  the  integrands  in  equation  (89) 
are  functions  of  z,  the  coordinate  along  the  propagation 
path.  Therefore,  in  order  to  determine  a  distortion  param¬ 
eter  for  thermal  blooming,  the  z  dependence  for  these  ele¬ 
ments  must  be  stated.  Of  course,  in  the  case  of  this 
study,  the  z  dependence  is  primarily  an  altitude  dependence. 

The  radius  of  the  beam  along  the  path  can  be  determined 
through  equation  (53)  after  considering  the  effects  of 
diffraction,  focusing,  jitter,  and  turbulence.  In  addition, 
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the  use  of  LOWTRAN  and  FASCODE  in  this  study  provides  values 
of  absorption  and  scattering  coefficients  for  molecular  and 
scattering  processes  at  various  layers  of  the  atmosphere. 
Thus,  total  attenuation  coefficients  and  total  absorption 
coefficients  can  likewise  be  determined  at  each  layer 
through  addition.  For  intermediate  values,  exponential 
interpolation  can  be  used.  Also,  LOWTRAN  will  provide 
values  of  pressure,  temperature  and  density  for  each  layer 
boundary  and  average  values  for  each  layer  in  the  assumed 
model  atmosphere.  Again,  interpolation  is  used  for  inter¬ 
mediate  altitudes.  Using  values  for  pressure , temperature , 
and  the  propagation  wavelength,  the  value  of  the  refractive 
index,  n(Z),  can  be  computed  through  equation  (68)  and  the 
differential  of  the  refractive  index  with  respect  to  temper¬ 
ature,  dn/dT,  can  likewise  be  computed  through  equation  (69) 
Both  equations  were  discussed  previously  under  turbulence 
effects . 

Expressions  for  the  specific  heat  at  constant  pressure, 

C  ,  and  the  nominal  crosswind,  V  ,  were  borrowed  from  Peck- 
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ham  and  Davis  (18:58).  We  will  assume  in  this  study  that 
a  horizontal  wind  always  exists  and  this  wind  continuously 
increases  with  altitude  until  the  tropopause,  at  which 
point,  the  average  wind  begins  to  decrease.  Peckham  and 
Davis'  curve  for  nominal  wind  velocity  versus  altitude  was 
based  on  data  contained  within  Reference  34.  A  plot  of  this 


curve  is  found  in  Figure  20.  Their  formulation  for  the 
specific  heat  curve  was  a  2nd  degree  polynomial  curve  fit 
with  temperature  based  on  data  in  Reference  35.  A  plot  of 
this  curve  is  found  in  Figure  21. 

Model  Resolution 

One  of  the  important  issues  in  the  development  of  the 
model  was  the  choice  of  resolution,  or  the  separation  be¬ 
tween  successive  wavelengths  which  would  be  compared.  The 
model  itself  will  pose  no  restrictions;  however,  a  trade¬ 
off  between  computer  time  and  the  discrimination  in  the 
output  data  was  necessary.  The  important  factors  bearing 
on  this  decision  will  be  discussed. 

As  described  earlier,  the  structure  of  the  molecular 
absorption  curve  is  very  complex  in  wavelength  regions 
longer  than  about  .55  micrometers.  Each  absorption  line 
has  a  finite  width.  At  sea  level,  collision  broadening 
mechanisms  dominate,  widening  the  average  absorption  line 
to  approximately  0.1  wavenumbers  cm  ^  (24:10).  At  higher 
altitudes,  the  linewidth  is  determined  primarily  by  the 
weaker  Doppler  broadening  to  be  approximately  0.01  cm  ^ 
(24:10).  FASCODE  samples  across  a  spectrum  at  a  specified 
interval.  Conceivably,  if  that  interval  is  greater  than 
the  average  linewidth  at  a  given  altitude,  the  entire  line 
contribution  would  be  missed.  If  fact,  the  original  FASCODE 


program  has  been  structured  to  sample  at  intervals  equal 
to  one  quarter  of  the  average  linewidth  for  this  reason 
(24:14). 

Also  to  be  considered  are  the  limitations  in  the  free 
electron  laser  device.  First,  the  relationship  between  the 
linewidth  of  the  laser  and  the  absorption  linewidth  must  be 
explored.  As  with  conventional  lasers,  the  free  electron 
laser  has  a  small  signal  gain  curve,  broadened  about  some 
central  value.  The  width  of  that  curve  is  based,  however, 
on  relativistic  dynamics  (11:11).  The  definition  of  which 
electrons  are  or  are  not  in  phase  can  only  be  determined 
within  the  length,  Lw,  of  the  wiggler.  Based  on  the  Un¬ 
certainty  Principle,  it  can  be  shown  that  the  fractional 
linewidth,  7v/v  ,  is  (11:11): 
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Typical  values  of  Lw=30 -meters  and  ^w=0.02  meters 
gives  a  linewidth  of  approximately  6.67  cm”'*'  for  visible 
frequencies.  However,  this  gain  curve  can  be  treated  as 
a  homogeneously  broadened  gain  curve  of  a  conventional 
laser;  the  laser  will  saturate  on  the  highest  gain  cavity 
mode  closest  to  the  center  of  the  gain  profile.  The 
emitted  linewidth  will,  therefore,  be  determined  by  the 
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cavity  mirror  stability  and  could  be  on  the  order  of  10"^ 
cm  \  many  orders  of  magnitude  narrower  than  the  average 
absorption  linewidth  in  the  atmosphere  (11:12).  Thus,  the 
laser  linewidth  does  not  restrict  the  model  resolution. 

Presumably,  with  a  little  judicious  cavity  adjust¬ 
ments,  any  wavelength  of  interest  could  be  selected.  How¬ 
ever,  an  open  question  remains  concerning  the  stability  of 
the  laser  wavelength  based  on  the  stability  of  the  electron 
beam.  No  reference  could  be  found  in  the  open  literature 
on  this  subject.  According  to  the  Free  Electron  Laser 
Program  office  at  Los  Alamos,  the  energy  stability  of  the 
electron  beam  can  be  controlled  by  feedback  devices  to  one 
part  in  10,000  or  better  (36).  As  a  gross  estimate  of  wave 
length  stability,  one  may  begin  with  equation  (10) 


Differentiating  with  respect  to  Y  ,  we  get, 

d A  .  =  Aw  d  y  (92) 

~T~ 

Assuming  Aw  =  0.02  meters,  *  ^  =  .  5  microns,  and  dy  =.0C005 
gives  a  spread  in  the  laser  output  of  approximately  2.0  cm" 
A  preliminary  test  run  of  the  model  indicated  that 
approximately  0.3  seconds  of  central  processing  time  was 
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required  for  each  wavelength  calculation.  Therefore,  a 
trade-off  between  capturing  the  majority  of  the  structure 
in  the  molecular  absorption  curve  and  computer  time  led  to 
a  decision  to  sample  in  intervals  of  0.2  wavenumbers.  It 
was  felt  that  the  structure  missed  by  this  choice  repre¬ 
sented  secondary  consequences.  Finally,  in  the  wavelength 
region  shorter  than  .55  micrometers  where  molecular  absorp 
tion  is  non-existent,  the  spectrum  will  be  sampled  at  a 
conservative  1.0  wavenumbers. 


Electromagnetic  radiation  incident  on  a  metallic  target 
will  either  be  absorbed  or  reflected.  Of  course,  for  laser 
weapon  system  effectiveness,  one  is  interested  in  the  degree 
of  absorption  of  laser  energy  into  the  target  material. 

The  free  electron  laser  normally  is  operated  to  produce 
a  long  pulsed  output  of  high  average  power  (2)  .  Therefore, 
the  primary  kill  mechanism  is  a  thermal  burn  through  of  the 
target's  outer  skin.  The  laser  energy  absorbed  in  the  skin 
is  converted  to  heat,  which  is  gradually  transported  through 
the  target  by  simple  conduction.  Such  thermal  kill  mechan¬ 
isms  can  be  assumed  to  be  linear  processes,  and  therefore, 
scale  directly  with  the  absorptivity  of  the  target.  In 
certain  regions  of  the  spectrum,  absorptivity  can  be  a  strong 
function  of  the  incident  laser  wavelength. 

Basic  Theory 

Consider  the  following  simple  modeL  A  monochromatic 
plane  wave  is  incident  perpendicular  to  a  smooth  conducting 
plane  surface.  The  electric  field  in  the  wave  causes 
oscillating  motion  in  the  loosely  bound  valence  electrons  of 
the  metal  at  the  same  frequency  of  the  incident  wave.  If 
totally  free  to  oscillate,  the  electrons  will  produce  a 
reflected  wave  of  the  same  magnitude,  exactly  180  degrees 
out  of  phase  with  the  incident  wave.  In  this  case  ,  no 
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radiation  is  absorbed  by  the  target  material  and  accounts  for 
the  highly  reflective  nature  of  metallic  conductors  like  sil¬ 
ver.  However,  viscous  damping  forces  arising  from  the  colli¬ 
sions  between  the  oscillating  electrons  and  the  atomic  lat¬ 
tice  of  the  metaL  will  restrict  the  movement  of  the  electrons, 
and  thus,  reduce  the  conductivity.  These  collisions  dissipate 
the  electrons'  energy  into  thermal  motion,  and  produce  heat 
within  the  target. 

The  above  discussion  constitutes  the  basic  assumptions 
behind  the  Drude  free  electron  model  (37:434).  The  theory 
pre-supposes  a  complex  index  of  refraction, 

n  ■  n  +  ik  (93) 

where  n  is  the  real  index  of  refraction  and  k  is  an  extinction 
coefficient  of  an  incident  wave.  The  extinction  component 
can  be  more  easily  seen  when  the  complex  index  of  refraction 
is  used  in  the  formula  for  the  wave's  electric  field. 

“koj  z  into 

E(z)  -  E(z  -  0)  e  e  ~C~  (94) 

The  Fresnel  expression  for  the  fraction  of  reflected 
radiation  from  a  homogenenous  surface  due  to  a  normally  in¬ 
cident  wave  is  given  by  the  following  (37:395). 
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For  optical  frequencies  metals  are  opaque;  therefore,  by 
Kirchoff's  Law  and  energy  observation  (37:395),  the  fraction 
of  radiation  reflected  plus  the  fraction  absorbed  must  equal 
one.  This  gives  an  expression  for  the  absorptance  of  normally 
incident  waves. 
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If  one  assumes  the  frequency  of  the  radiation  is  less 
than  the  frequency  of  the  induced  electron  collisions  within 
the  target  material,  then  the  currents  in  the  metal  are  in 
phase  with  the  electric  field  of  the  radiation.  In  this  case, 
the  absorption  depends  on  just  the  incident  wavelength,  X,  and 
the  conductivity  of  the  material,  o.  This  expression  is  wide¬ 
ly  known  as  the  Hagen-Rubens  formula  (37:437). 


(97) 


where  c  is  the  speed  of  light,  and  is  the  permittivity  of 
free  space. 

By  this  expression,  the  absorptivity  of  a  metallic  tar¬ 
get  scales  as  A"^  and  therefore,  tends  to  favor  shorter 
wavelengths  for  a  laser  weapon  system  analysis.  For  pure 
metals  at  room  temperature,  equation  (97)  is  usually  not 
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valid  in  the  visible  portion  of  the  spectrum.  Since  visible 
frequencies  are  generally  greater  than  the  induced  electron 
collision  frequencies,  this  model  would  tend  to  over  estimate 
the  absorptance  of  the  metal.  However,  for  alloys  of  common 
aerospace  materials  where  impurities  exist  in  the  atomic 
lattice,  the  experimental  data  suggests  that  this  expression 
may  be  quite  appropriate. 

Figure  22  shows  the  spectral  absorptivities  for  four 
common  alloys  used  in  satellites,  missiles,  and  aircraft  as 
a  function  of  wavenumber.  In  the  visible  region,  10000  to 
29000  cm"^,  the  absorptivity  does  seem  to  indeed  scale  as 
equation  (97).  This  assumption  will  be  tested  later  in 
this  chapter. 

However,  this  simple  scaling  is  far  from  explaining  the 
coupling  of  laser  radiation  to  a  target.  The  spectral  ab¬ 
sorption  curve  is  greatly  modified  by  the  existance  of  re¬ 
sonance  wavelengths.  This  phenomenon  is  not  unlike  the 
natural  frequencies  of  an  oscillator  system.  This  behavior 
can  be  observed  in  the  Aluminum-2024  curve  in  figure  22 
where  an  absorption  peak  occurs  in  the  visible  spectrum  at 
approximately  12000  cm”^(.83  micrometers). 

In  addition,  the  surface  conditions  of  metallic  targets 
play  a  dominant  role  in  their  radiative  properties.  The 
degree  of  surface  roughness  can  be  severe  enough  to  require 
a  statistical  treatment  in  order  to  explain  the  optical  be¬ 
havior  of  the  surface.  If  the  surface  irregularities  are 
much  smaller  than  the  wavelength  of  light,  diffraction 
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theory  can  be  used.  In  this  case, it  has  been  shown  that  the 
surface  reflectance  becomes  an  exponential  function  of  the 
square  of  the  laser  frequency  (39:18).  Another  very  import¬ 
ant  surface  effect  is  the  existence  of  thin  films,  either  in¬ 
troduced  by  man  or  the  result  of  oxide  growths.  These  films 
act  much  like  interference  filters  where  certain  wavelengths 
are  selectively  reflected  depending  on  the  thickness  of  the 
film  (39:18).  This  effect  is  also,  of  course,  very  dependent 
on  the  angle  of  incidence  of  the  laser  wave.  Finally,  the 
spectral  absorptance  curve  of  a  meterial  is  often  dominated 
by  other  factors  such  as  the  thermal  history  of  the  material, 
chemical  preparations,  alloying,  lattice  orientation,  im¬ 
purities,  and  the  presence  of  paint  (39:18). 

Being  confronted  with  these  complexities  proved  to  be  a 
difficult  problem  for  this  study,  since  it  was  desired  not 
to  make  some  specific  assumptions  concerning  the  target  com¬ 
position.  The  use  of  actual  target  spectral  data  would  quick 
ly  classify  this  study  as  well.  All  this  led  to  the  decision 
to  scaLe  the  wavelength  dependence  of  the  target  interaction 
with  the  simplified  Crude  free  electron  theory  as  in  equation 
(97).  For  this  study,  we  will  define  a  value  of  target 
'efficiency'  between  zero  and  one.  The  best  efficiency  will 
occur  at  the  short  wavelength  end  of  the  spectrum  of  interest 
for  this  study,  i.e.  29000  cm”^  or  .3448  micrometers.  For 
longer  wavelengths,  efficiency  values  will  be  scaled  pro¬ 
portionally  downward  from  one.  The  expression  used  in  the 


program  will  be  the  following. 


Target  efficiency  * 


.3448  (microns) 
X  (microns ) 


% 


(98) 


The  program  will  be  structured  in  such  a  way  that  more  real¬ 
istic  target  data  can  easily  be  inserted  at  a  later  time. 


Discussion  of  the  Approach 

The  approach  we  have  selected  is  based  on  two  rather 
questionable  assumptions:  1)  the  effectiveness  of  a  laser 
weapon  is  more  a  function  of  the  bulk  properties  of  the  tar¬ 
get  material  rather  than  the  surface  properties  of  the  tar¬ 
get  and  2)  the  bulk  properties  can  be  approximated  through 
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the  Drude  free  electron  model  and  a  X  scaling.  In  this 
section,  we  will  further  examine  the  validity  of  these  as¬ 
sumptions. 

An  Air  Force  Weapons  Laboratory  technical  report  de¬ 
scribes  two  distinct  mechanisms  for  the  removal  of  paint  from 

a  metallic  target  surface  (40:3).  At  relatively  low  laser 

2 

intensities  on  the  order  of  kilowatts/cm  ,  the  paint  appar¬ 
ently  decomposes  and  chars,  leaving  the  metal  target  blackened, 
and  thus,  highly  absorbing  (40:25).  In  this  case,  the  exis- 
ence  of  a  paint  actually  enhances  target  coupling.  However, 
at  higher  laser  intensities  (approximately  10  KW/cm  ),  the 
paint  becomes  quickly  vaporized  and  leaves  behind  little  res- 


idue.  In  this  case,  the  average  absorptance  of  the  painted 
target  approaches  the  absorptance  of  the  unpainted  metal 
(40:26). 

The  second  assumption  we  have  made  concerning  the  scal¬ 
ing  of  target  interactions  can  be  justified  based  on  experi¬ 
mental  data.  Such  data  on  the  radiative  properties  of  mate¬ 
rials  can  be  found  in  Reference  39  for  common  aerospace  al¬ 
loys.  We  selected  two  common  alloys:  Aluminum- 2024 ,  a 
wrought  alloy  with  copper  as  the  principal  alloying  agent 
(39:27),  and  Auminum-7075 ,  also  a  wrought  alloy  with  zinc  as 
the  alloying  agent  (39:114).  Each  experiment  cataloged  con¬ 
sisted  of  absorption  measurements  at  various  wavelengths  for 
the  given  sample  of  material.  These  experiments  represented 
a  wide  variety  of  surface  preparations  and  thermal  histories 
for  the  material  samples  as  well.  The  raw  data  from  these 
experiments  are  found  on  pages  37,  42,  109,  115,  119,  and 
136  of  Reference  39.  For  our  analysis,  the  raw  absorption 
coefficients  resulting  from  each  experiment  were  normalized 
to  one  at  the  value  corresponding  to  .35  micrometers.  In 
all,  49  data  points  were  used. 

In  this  analysis,  we  assumed  a  statistical  model  where 
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the  normalized  absorption  data  is  proportional  to  A  -.  The 
data  was  transformed  in  order  to  achieve  a  linear  relationship 
between  the  independent  variable  of  wavelength  and  the,  now 
transformed,  dependent  variable  of  absorption.  The  data  was 


then  inserted  into  a  simple  univariate  linear  regression  rou¬ 
tine  . 

The  two  variables  were  found  to  be  well  correlated  under 
this  assumed  model.  The  correlation  coefficient  from  the  re¬ 
gression  was  .794,  indicating  that  79.4%  of  the  variation  of 
the  data  can  be  explained  through  our  assumed  model.  Much 
of  the  variation  in  the  data  is  from  the  resonance  absorption 
feature  in  aluminum  at  approximately  12000  cm”  (.83  microm¬ 
eters).  The  linear  regression  line  between  the  transformed 
variables  was  found  to  be: 

Absorption  *  -.0608  +  3.1604  X  (micrometers)  (99) 

This  relationship,  transformed  back  to  the  original  variables, 
is  plotted  in  Figure  23  along  with  the  experimental  data  from 
which  it  was  determined.  The  data  is  displayed  in  wavenumbers 
to  keep  consistency  with  other  spectral  plots  in  this  report. 
The  maximum  deviation  between  the  regression  curve  and  our 
target  efficiency  expression,  equation  (98),  is  3.6%. 

This  analysis  is  far  from  conclusively  showing  that  all 
laser  beam/target  interactions  can  be  described  through  equa¬ 
tions  (97).  The  subject  of  target  energy  coupling  is  extreme¬ 
ly  complex,  and  perhaps,  can  only  be  fully  treated  empirically 
for  a  specific  target  material. 


The  theory  developed  in  the  previous  chapter  provided  the 
foundation  for  the  model.  This  chapter  outlines  the  actual 
program  which  was  used  to  examine  the  relationship  between 
wavelength  and  the  ground  based  laser  weapon  system  efficiency 
to  include  propagation,  device,  and  target  coupling  efficien¬ 
cies.  To  those  readers  not  interested  in  the  actual  program, 
this  chapter  may  be  skipped  without  loss  of  continuity. 

The  program  was  coded  in  FORTRAN  V  for  use  on  a  Control 
Data  Corporation  Cyber  6600  series  computer  with  the  NOS/BE  1 
operating  system.  The  overall  program  flow  requires  the  execu 
tion  of  five  separate  programs  as  depicted  in  Figure  24.  At 
Wright-Patterson  AFB,  all  of  the  molecular  absorption  line 
data  for  the  program  FASCODE  is  contained  on  magnetic  tape. 
This  is  preferred  for  two  reasons:  1)  to  reduce  storage  re¬ 
quirements  on  the  computer,  and  2)  to  minimize  input/output 
operations  during  execution  of  FASCODE.  Thus,  the  first  step 
is  to  read  this  data  from  tape  to  disk  computer  memory.  This 
data  is  independent  of  the  atmospheric  model  or  transmission 
geometry  chosen  by  the  user;  therefore,  this  data  can  be 
transferred  once  for  a  portion  of  the  spectrum  and  accessed 
during  multiple  runs  of  FASCODE.  The  Air  Force  Geophysics  Lab 
oratory  makes  available  programs  for  this  purpose,  thus  assur¬ 
ing  compatibility  with  the  main  program,  FASCODE  (24:68). 


The  program  FASCODE  was  modified  to  provide  molecular  ab¬ 
sorption  optical  depths  for  each  of  32  layers  of  the  atmosphere 
so  that  thermal  blooming  effects  could  be  calculated.  The  mod¬ 
ifications  to  FASCODE  are  discussed  in  Appendix  B.  The  term 
optical  depth  is  that  value  where  the  laser  beam  is  attenuated 
by  a  factor  of  1/e.  It  is  computed  by  simply  multiplying  the 
attenuation  coefficient,  a,  described  by  equation  (63),  by  the 
path  length,  Az.  In  this  case,  the  path  length  is  the  width 
of  each  altitude  layer. 

Unfortunately  for  this  study,  FASCODE  computes  optical 
depths  for  all  desired  wavenumbers  for  an  entire  layer  at  a 
time.  However,  our  main  program,  WAVLEN,  requires  absorption 
data  for  all  layers,  one  wavenumber  at  a  time.  This  incompat¬ 
ibility  required  that  a  program  to  sort  the  FASCODE  absorption 
data  to  a  usable  format  also  be  developed.  This  program,  SORT, 
is  included  in  Appendix  C.  The  sorting  process  produces  a  data 
file  of  optical  depths  called  FCDATA. 

The  program  LOWTRAN  calculates  its  optical  depths  in  a 
format  compatable  with  the  main  program  and  thus  requires  no 
additional  data  manipulation.  As  modified,  LOWTRAN  creates  a 
data  file  called  LTDATA  which  contains  the  attenuation  coef¬ 
ficients  of  Rayleigh  scattering,  aerosol  absorption,  aerosol 
scattering,  continuum  absorption,  and  ozone  molecular  absorp¬ 
tion  for  each  altitude  layer.  The  modification  of  LOWTRAN  is 
described  in  Appendix  A. 

Both  FASCODE  and  LOWTRAN  use  identical  atmospheric  models 


in  their  calculations.  These  values  of  temperature,  pressure, 
and  density  at  the  various  altitudes  are  also  required  in  the 
main  program,  WAVLEN.  It  was  arbitrarily  decided  that  LOWTRAN 
would  provide  this  data  by  creating  a  data  file  called  MODEL 
during  the  execution  of  LOWTRAN. 

These  three  data  files  provide  input  to  the  main  program, 
WAVLEN.  Since  this  program  is  the  guts  of  the  computer  model, 
it  will  be  discussed  in  some  detail.  A  listing  of  program 
WAVLEN  is  found  in  Appendix  D. 

The  purpose  of  the  program  is  to  calculate  a  weapon  system 
efficiency  for  a  range  of  wavenumbers  separated  by  a  specified 
interval.  Through  inspection  of  the  program  output,  an  optim¬ 
um  efficiency  is  then  determined.  Therefore,  the  structure  of 
the  program  is  identified  by  one  large  wavenumber  loop.  A  flow 
chart  for  program  WAVLEN  is  depicted  in  Figure  25  for  refer¬ 
ence.  The  program  requires  ten  input  parameters  for  operation. 
These  are: 

AVGPWR  --  Average  laser  output  power  in  Watts 

RO  --  The  radius  of  the  output  optics  in  meters 
VI  --  Beginning  wavenumber  to  be  calculated  in  cm-1 

V2  —  Ending  wavenumber  to  be  calculated  in  cm-1 

RANGE  --  Range  from  the  laser  to  the  target  in  kilometers 

FOCRNG  —  Focal  range  of  the  output  optics  in  kilometers 

BEAMQ  --  Beam  quality,  times  diffraction  limited 
JITTER  —  Mean  single  axis  jitter  in  microradians 
HOL  --  Height  above  mean  sea  level  of  laser  in  kilometers 
COMPEN  --  Percentage  of  the  turbulence  induced  beam 
radius  corrected  by  adaptive  optics 

All  distance  measurements  in  the  program  are  converted  to 
meters  for  consistency.  In  fact,  MRS  units  are  used  exten- 


sively.  One  of  the  first  functions  performed  by  WAVLEN  is  to 
read  in  the  atmospheric  data  from  the  data  file  MODEL.  Sub¬ 
routine  MDLATM  performs  this  operation.  Values  of  pressure 
in  millibars,  temperature  in  degrees  Kelvin,  and  density  in 
grams/cubic  cm  for  both  the  altitude  boundaries  and  average 
values  for  each  altitude  layer  are  read.  Density  values  are 
converted  to  kilograms/cubic  meters,  as  well. 

The  next  operation  begins  the  calculation  of  Yura's  co¬ 
herence  diameter,  Pq.  This  is  done  now,  outside  of  the  main 
wavelength  loop,  to  save  computer  time  since  much  of  the  cal¬ 
culations  are  wavelength  independent.  Two  different  models 
are  used  for  the  refractive  index  structure  parameter  varia¬ 
tion  with  altitude.  From  ground  level  to  1000  meters,  corre¬ 
sponding  to  the  first  temperature  inversion,  Wyngaard's  model 
is  used.  The  environmental  conditions  chosen  for  this  model 
represent  mid-morning  or  mid-afternoon  sun  angles  during  a 
clear  day.  Further,  the  height  of  the  laser  is  assumed  to  be 
50  meters  above  dry,  open  fields.  The  summation  for  Yura's 
coherence  diameter  proceeds  in  50  meter  increments  for  these 
first  1000  meters.  For  this  model,  a  value  for  the  change  in 
refractive  index  for  a  change  in  temperature  is  required  and 
is  provided  by  the  external  function,  DNDT.  The  arguments  for 
this  function  are  local  temperature,  pressure,  and  wavelength. 
An  average  value  of  .6  microns  is  used.  Intermediate  values 
of  temperature  are  found  through  linear  interpolation.  How- 
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Figure  25.  Flow  Chart  of  Program  WAVLEN 
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ever,  exponential  interpolation  is  used  to  determine  values  of 
pressure.  Here,  we  have  assumed 


|  «  DENSITY  «  exp  (-z/SCALE  HEIGHT)  (100) 

For  any  generic  exponential  function,  f(z),  the  graph  of 
the  natural  logarithm  of  f(z)  is  a  straight  line.  Therefore, 
intermediate  values  between  z*a  and  z*b  can  be  determined  by 
first  employing  linear  interpolation. 

In  Hz)  -  In  f (a)  ♦  (§-=-f  j  ln  ($£})  <101> 

Taking  the  exponential  of  both  sides  of  equation  (101)  yields: 


Above  1000  meters  to  a  height  of  25  kilometers,  Hufnagel's 
model  of  the  refractive  index  structure  parameter  is  employed. 
Here,  the  summation  to  determine  Pq  is  in  increments  of  1000 
meters . 

The  program  WAVLEN  next  enters  the  wavenumber  loop  and 
calls  subroutine  FCREAD  to  read  in  the  molecular  absorption 
information  from  the  data  file  FCDATA.  The  wavenumber  read 
by  FCREAD  determines  that  wavenumber  value  to  be  used  through¬ 
out  the  rest  of  the  program.  A  provision  is  made  in  the  program 
to  accomodate  runs  which  transition  over  the  end  of  the  spectral 


region  covered  by  FASCODE  (approximately  18000  cm-^).  When 
this  point  is  reached,  the  variable  IFLAG  is  set  to  one  and 
the  program  continues  with  a  specified  wavenumber  interval 
(we  have  chosen  1.0  cm-1)  and  using  just  the  LOWTRAN  data  for 
input. 

The  program  passes  the  wavenumber  read  by  FCREAD  when  it 
next  calls  subroutine  LTREAD.  This  subroutine  reads  in  the 
information  contained  in  the  data  file  LTDATA.  This  data  is 
separated  by  5.00  cm”'*'  wavenumber  intervals.  With  one  excep¬ 
tion,  simple  linear  interpolation  is  used  to  find  the  layer  at¬ 
tenuation  for  the  wavenumber  of  interest.  In  the  case  of  Ray¬ 
leigh  scattering,  however,  the  interpolation  is  scaled  by  wave- 
number  values  to  the  fourth  power.  If  an  end-of-file  condition 
is  experienced  on  the  file  LTDATA,  the  program  sets  IFLAG  to 
two  and  processing  ceases. 

To  save  computer  time,  WAVLEN  makes  a  quick  determination 
of  atmospheric  transmission  and  thermal  blooming  to  see  if  it 
is  worthwhile  proceeding  with  further  calculations  for  the  cur¬ 
rent  wavenumber.  The  total  transmission  is  determined  by  ad¬ 
ding  the  contributions  from  all  of  the  atmospheric  attenuation 
processes.  A  rough  value  of  Gebhardt's  thermal  blooming  dis¬ 
tortion  parameter,  based  on  conditions  in  the  first  altitude 
layer,  is  also  calculated.  If  either  the  transmission  is  less 
than  10%  or  the  thermal  blooming  distrotion  parameter  is  great¬ 
er  than  30  (corresponding  to  a  98%  irradiance  loss  in  the  laser 
beam),  then  the  overall  efficiency  of  the  weapon  system  is  set 
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to  zero. 


Following  completion  of  the  wavelength  dependent  portion 
of  the  computation  of  Yura's  coherence  diameter,  Pq,  the  pro¬ 
gram  begins  the  thermal  blooming  calculations.  The  double 
integration  involved  in  computing  the  thermal  blooming  dis¬ 
tortion  parameter,  NI,  requires  two  altitude  intervals  over 
which  to  sum.  In  the  outer  summation,  the  interval  is  simply 
the  width  of  the  individual  altitude  layers.  For  the  inner 
summation,  however,  each  interval  is  reduced  by  a  factor  of 
one  fourth.  Tests  showed  that  there  is  no  problem  of  con¬ 
vergence  of  the  integral  with  these  chosen  intervals;  reduc¬ 
ing  the  intervals  did  not  alter  the  final  result. 

As  with  the  previous  calculation  of  Yura's  coherence 
diameter,  intermediate  values  of  temperature  are  found  through 
linear  interpolation  and  pressure  through  exponential  interpo¬ 
lation.  Once  these  local  values  of  temperature  and  pressure 
are  determined,  they  are  then  used  to  scale  the  density. 

The  determination  of  the  absorption  and  attenuation  coef¬ 
ficients  for  the  thermal  blooming  calculations  requires  some 
explanation.  The  optical  depths  computed  by  FASCODE  and  LOWTRAN 
both  represent  mean  values  for  a  given  altitude  layer.  Since 
the  optical  depths  are  proportional  to  the  number  density  of 
absorbers  or  scatterers,  which  in  turn  generally  scale  expon¬ 
entially  with  altitude,  the  mean  value  of  optical  depth  does 
not  coincide  with  the  middle  of  the  altitude  layer.  If  we  as- 
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sume  a  constant  scale  height  for  the  atmosphere  of  7.0  kilom¬ 
eters  and  borrow  the  concept  of  'expected'  or  mean  value  from 
statistics  (41:135),  we  see  that 

1 

f  -Z/7.0 

Z  -  I  Ze  dZ  =  .46  (103) 

0 

Therefore,  the  mean  value  of  the  optical  depth  occurs,  not  at 
the  midpoint,  but  rather  at  an  altitude  which  is  46%  from  the 
bottom  of  the  layer.  Using  these  points  as  references,  we 
can  nov  determine  intermediate  values  of  optical  depth  with 
exponential  interpolation.  Finally,  to  convert  the  optical 
depth  to  an  absorption  or  attenuation  coefficient,  one  must 
divide  by  the  width  of  the  altitude  layer. 

To  complete  the  thermal  blooming  calculations,  function 
REF1DX  provides  the  local  index  of  refraction  given  the  local 
pressure  and  temperature  and  the  propagating  wavelength  and 
function  VX  supplies  a  nominal  value  of  the  horizontal  wind 
for  a  given  altitude.  Once  the  thermal  blooming  distortion 
parameter,  NI,  is  computed,  it  becomes  the  argument  to  func¬ 
tion  IREL  to  determine  the  irradiance  reduction  factor  for 
the  beam. 

2 

The  1/e  radius  of  the  laser  beam  at  any  point  along  the 
propagation  path  is  determined  by  subroutine  RADIUS.  It  com¬ 
bines,  in  a  root-sum-squared  technique,  the  characteristic 


radii  resulting  from  the  effects  of  turbulence,  jitter,  and 
diffraction  and  focusing.  For  focused  laser  beam,  the  actual 
focal  range  of  the  output  optics  is  used;  however,  if  a  colli¬ 
mated  laser  beam  is  desired,  the  user  should  input  a  value  of 

zero  for  the  focal  range.  The  program  then  sets  the  focal 
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range  to  infinity  (actually  10  meters)  for  all  beam  radius 
calculations.  For  turbulence,  the  program  includes  the  ability 
to  specify  the  percentage  that  atmospheric  compensation  by  a- 
daptive  optics  systems  reduces  the  turbulence  induced  beam 
radius . 

The  last  portion  of  program  WAVLEN  determines  the  different 
efficiencies  at  a  given  wavenumber  and  outputs  the  results.  A 
propagation  efficiency  is  computed  according  to  equation  (54). 
The  free  electron  laser  efficiency  with  energy  recovery  is  sup¬ 
plied  by  subroutine  FELOPT  and  a  target  coupling  efficiency  is 
provided  by  function  TGTEFF.  These  three  efficiencies  are  then 
multiplied  together  to  arrive  at  an  overall  weapon  system  ef¬ 
ficiency  for  that  wavenumber.  The  program  automatically  keeps 
track  of  the  best  propagation  efficiency  and  best  overall  ef¬ 
ficiency  over  the  run. 

Program  Limitations 

Aside  from  the  basic  assumptions  in  the  theory  and  the 


assumptions  we  have  chosen  to  make,  the  program  possess  certain 
limitations.  Even  though  LOWTRAN  and  FASCODE  can  compute  op¬ 
tical  depths  for  slant  paths  through  the  atmosphere,  WAVLEN 


can  only  accomodate  vertical  propagation  paths.  Primarily, 
this  is  due  to  the  many  altitude  dependent  parameters  used  in 
the  turbulence  and  thermal  blooming  calculations.  Also,  the 
program  is  structured  to  use  the  entire  atmosphere,  i.e.  the 
selected  range  to  the  target  must  be  greater  than  100  kilometers. 
Finally,  the  height  of  the  laser  must  correspond  to  one  of  the 
boundaries  of  the  altitude  layers.  In  essence,  this  requires 
that  the  laser  height  be  specified  to  the  nearest  kilometer 
above  sea  level.  The  parameter,  COMPEN,  which  adjusts  the  ef¬ 
fects  of  atmospheric  turbulence  on  the  beam  radius,  is  strictly 
a  scaling  factor  and  should  not  be  confused  with  degrees  of  com¬ 
pensation  quoted  for  actual  adaptive  optics  systems. 

It  should  be  emphasized  that  the  efficiency  values  derived 
by  the  program  do  not  represent  actual  values.  In  some  cases, 
proportionality  constants  have  been  excluded.  The  purpose  of 
the  study  is  to  assess  the  relative  efficiencies  between  wave¬ 
lengths  and  not  to  determine  the  actual  values. 
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VI  .  Analysis  Results 

The  model  developed  for  this  study  and  presented  in  the 
preceeding  chapters,  was  exercised  in  order  to  examine  the 
interactions  between  the  wavelength  dependent  relationships 
in  the  system  and  to  select  a  'best'  wavelength  for  operation 
This  chapter  presents  the  results  of  this  analysis. 

Selection  of  Parameters  for  the  Baseline  Case 

A  baseline  case  was  used  as  a  point  of  departure  for  this 
analysis  and  to  validate  some  of  the  factors  within  the  model. 
Therefore,  the  selection  of  input  parameters  for  this  base¬ 
line  case  was  somewhat  critical.  This  section  will  discuss 
that  selection  process. 

A  mid-latitude,  summer  atmospheric  model  was  used  for 
the  baseline  calculations.  In  addition,  a  Rural  boundary 
layer  aerosol  model  with  a  23  km  visibility  was  initially 
assumed. 

An  important  consideration  for  this  analysis  is,  of 
course,  the  wavelength  region  to  be  examined.  This  region 
must  be  within  the  expected  capabilities  of  the  free  electron 
laser  device  and  should  include  areas  where  the  atmospheric 
attenuation  is  low.  Presumably,  this  region  falls  in  or 
near  the  visible  portion  of  the  spectrum.  We  have  chosen 
the  wavelength  range  from  1.14  micrometers  to  .34  micrometers 
(8750.0  cm”l  to  29000.0  cm”^  for  this  study.  At  wavelengths 


shorter  than  .34  micrometers , it  is  expected  that  Rayleigh 
scattering  and  the  onset  of  significant  ozone  attenuation  in 


the  ultraviolet  region  will  combine  to  negate  any  advantage 
to  shorter  wavelengths.  Likewise,  at  wavelengths  longer  than 
1.14  micrometers,  the  spectral  atmospheric  transmission  curve 
enters  a  strongly  absorbing  region  due  to  water  vapor  (see 
Figure  12)  and  would  preclude  laser  operation  there. 

The  selection  of  laser  power  and  radius  of  the  output 
optics  went  hand  in  hand.  In  advocating  short  wavelengths 
laser  systems  for  the  military,  Canavan  uses  the  scaling 
parameter  of  'brightness'  to  assess  the  capabilities  of 
different  laser  systems  (42:3).  Brightness  is  the  power 
per  solid  angle  illuminated  by  the  laser  system  and  can  be 
approximated  by  the  following  expression. 

2 

Brightness  =  ^2 -  (104) 

Az 

where  P  is  the  laser  power,  D  is  the  diameter  of  the  output 

optics,  and  A  is  the  laser  wavelength.  Ignoring  atmospheric 

losses,  Canavan  claims  that  in  order  to  accomplish  a  strategic 

defense  mission  against  ballistic  missiles,  a  laser  system 

21 

would  have  to  produce  a  brightness  on  the  order  of  10 
watts/steradian  (42:7).  For  wavelengths  in  the  visible,  this 
brightness  can  be  achieved  through  a  10  megawatt  laser  with 
5  meter  optics. 

Such  a  combination  of  power  and  optics  could  easily 
accomplish  an  anti-satellite  mission  for  any  conceivable 
orbit  altitude.  However,  the  ballistic  missile  defense 
mission  would  require  relay  mirrors  in  1000  to  2000  km  orbits 


(43:79).  If  one  restricts  oneself  to  shooting  through  the 
atmosphere  at  zenith  angles  no  greater  than  45  degrees,  simple 
geometry  tells  us  that  two  such  mirrors  would  be  required  and 
would  result  in  a  total  propagation  distance  from  laser  to 
target  on  the  order  of  20,000  km.  Therefore,  for  this  base¬ 
line  study,  a  selected  range  of  20,000  km  appeared  appropriate. 
In  addition,  a  focal  range  of  20,000  km  was  chosen,  equaling 
approximately  one  half  of  the  average  Rayleigh  range. 

The  choices  of  1.2  for  beam  quality,  0.2  microradians  for 
jitter,  1000  meters  for  the  mean  sea  level  altitude  of  the 
laser  and  .85  for  the  degree  of  turbulence  compensation  were, 
admittedly,  somewhat  arbitrary.  The  beam  quality  for  a  free 
electron  laser  is  expected  to  be  quite  good,  and  the  chosen 
value  represents  this  assumption.  The  free  electron  laser 
does  not  contain  large  mechanical  moving  parts,  and  therefore, 
is  not  expected  to  suffer  significantly  from  jitter.  The 
value  of  .2  microradians  represents  an  optimistic  jitter 
component.  Presumably,  the  laser  weapon  would  not  be  placed 
at  a  sea  level  altitude,  and  therefore,  1000  meters  might  re¬ 
present  a  suitable  starting  point.  As  mentioned  earlier, 
the  factor  used  in  this  model  for  turbulence  compensation  is 
not  to  be  compared  to  actual  compensation  factors  quoted  for 
adaptive  optics  systems;  it  is  simply  treated  as  a  scaling 
parameter.  In  the  visible  region,  approximately  85%  of  the 
total  turbulence  induced  beam  spread  is  due  to  the  instantan¬ 
eous  ,  small  scale  beam  spread  component  (44:30).  A  choice 
of  .85  for  the  baseline  case,  thus,  represents  outstanding 


compensation  for  this  small  scale  component  and  places  the 
turbulence  induced  beam  radius  on  the  same  order  of  magnit 
as  those  of  diffraction  and  jitter  effects. 


Therefore,  to  summarize,  the  baseline  case  for  this 


analysis  will  include  the  following. 


Average  laser  power 
Radius  ofoutput  optics 

Initial  wavenumber 

Final  wavenumber 
Propagation  range 
Focus  range 
Beam  quality 
Single  axis  jitter 
Height  of  laser  (MSL) 
Turbulence  compensation 


10  megawatts 
2.5  meters 

8750.0  cm"1 

29000.0  cm'1 
20,000  km 
20,000  km 
1.2 

.2  microradians 
1000  meters 
.85 


Baseline  results 

The  following  12  pages  (Figures  26  through  37)  contai 
the  plotted  results  of  the  baseline  case.  The  plots  have 
been  mistakenly  normalized  such  that  an  efficiency  value  o 
one  corresponds  to  18200.0  cm  1.  Although  this  results  in 
displayed  efficiencies  greater  than  one,  the  relative  effi 
ciencies  between  wavenumbers  remains  the  same. 

As  expected,  the  structure  of  the  efficiency  curve  is 
dominated  by  the  variations  caused  by  atmospheric  molecula 
absorption.  Most  of  the  absorption  lines  in  the  region  be 
tween  8750  and  13000  cm"1  are  caused  by  water  vapor  and 
carbon  dioxide,  and  in  some  places,  completely  negates 
laser  operation  there. 
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Figure  26.  Overall  System  Efficiency  (  8750  -  9500  cm 
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Figure  33.  Overall  System  Efficiency  (  15500  -  16500  cm 


WAVENUMBER 


The  optimum  wavenumber  for  this  baseline  scenario 

occurs  at  14347.2  cm  ^  (  .697  micrometers)  with  a  computed 

efficiency  of  .5079%.  Yet,  the  term  'optimum'  must  be  used 

with  caution.  Efficiencies  within  1%  of  this  optimum  value 

can  be  found  from  12476.0  to  15225.6  cm"^  (approximately 

.657  to  .801  micrometers).  Of  course,  the  cumulative  errors 

and  basic  assumptions  within  the  model  far  exceed  1%.  The 

relatively  flat  nature  of  the  top  line  structure  or  'envelope 

_  1 

of  the  efficiency  curve  between  10000  and  17500  cm  was  not 
expected.  It  would  appear  at  first  glance  that  the  compe¬ 
tition  among  the  wavelength  dependent  relationships  produces 
an  overall  system  efficiency  which  changes  very  little  over 
a  relatively  broad  region.  Any  wavelength  selected  within 
this  region,  as  long  as  it  avoided  a  molecular  absorption 
feature,  would  essentially  be  equally  efficient.  In  the 
next  sections,  we  will  proceed  to  examine  these  relation¬ 
ships  further  in  an  attempt  to  understand  exactly  what  is 
driving  the  top  level  structure  of  the  efficiency  curve. 

Thermal  Blooming.  With  the  parameters  chosen  for  the 
baseline  case,  thermal  blooming  effects  were  not  significant. 
Tests  runs  were  accomplished  where  the  laser  power  was  in¬ 
creased  in  10  megawatts  increments.  Not  until  a  laser  power 
of  50  megawatts  was  reached  did  the  results  change  apprecia¬ 
bly.  Even  then,  only  regions  of  significant  molecular 
absorption  were  affected.  In  no  instant  did  thermal  blooming 
affect  the  overall  efficiency  when  the  wavelength  did  not 


coincide  with  a  molecular  absorption  feature.  The  aerosol 
absorption  in  the  atmosphere  computed  by  LOWTRAN  proved  to 
be  a  very  small  percentage  of  total  atmospheric  absorption, 
Therefore,  aerosols  did  not  contribute  to  heating  the  at¬ 
mosphere  to  any  significant  degree.  Another  reason  for  the 
minor  impact  of  thermal  blooming,  even  within  absorbing  re¬ 
gions,  was  perhaps  the  model  we  selected  for  the  local  wind 
At  the  1000  meter  altitude  of  the  laser,  the  local  wind  was 
computed  to  be  11.2  meters/sec  or  25  miles  per  hour.  This 
value  seems  rather  large  and  could  perhaps,  in  and  of  itsel 
have  prevented  the  onset  of  thermal  blooming  in  those  cases 
involving  modest  laser  powers  and  weakly  absorbing  atmosphe 
ric  conditions. 

Atmospheric  Transmission.  The  individaul  optical  dept 
due  to  aerosols,  ozone  and  Rayleigh  scattering  were  extrac 
ted  during  program  execution  and  summed  over  the  applicable 
atmospheric  layers.  With  these  values,  a  transmission  fac¬ 
tor  attributable  to  each  effect  was  determined.  These  were 
used  to  examine  the  overall  structure  of  the  efficiency 
curve  and  to  validate  this  portion  of  the  model  as  well. 

The  individual  spectral  transmission  curves  for  these 
three  effects  are  found  in  Figures  38,  39,  and  40.  As  ex¬ 
pected,  aerccol  transmission  tended  to  favor  the  longer 
wavelengths  (lower  wavenumbers).  This  factor  proved  to  be 
a  rather  significant  issue,  varying  by  a  factor  of  two  thir 
over  the  spectral  range  of  the  study.  The  slight  undulatio 
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in  the  curve  are  presumably  caused  by  LOWTRAN's  use  of 
specific  aerosol  size  distributions  in  computing  attenuation 
coefficients,  as  explained  in  Chapter  III.  Interpolation 
within  both  LOWTRAN  and  our  program  undoubtedly  lead  to  fur¬ 
ther  inaccuracies. 

The  ozone  transmission  curve  reflects  the  weak  absorp¬ 
tion  feature  of  ozone  in  the  midportion  of  the  visible  spec¬ 
trum.  This  curve  is  identical  to  the  curve  in  Figure  16. 

This  absorbing  region  causes  the  shallow  yet  significant, 
dip  in  the  computed  efficiency  curves  between  15000  and 
18200  cm  ^ .  Also,  one  can  note  the  beginnings  of  the  ultra¬ 
violet  ozone  absorption  in  the  far  right  of  the  curve  in 
Figure  39.  Since  the  majority  of  ozone  attenuation  occurs 
at  higher  altitudes,  the  height  of  the  laser  is  expected  to 
have  little  bearing  on  this  curve. 

Finally,  the  Rayleigh  scattering  curve  depicts  the 
strong  wavelength  dependence  of  this  effect.  In  fact,  it 
can  be  shown  that  in  the  high  wavenumber  region  of  the  spec¬ 
trum,  Rayleigh  scattering  is  indeed  the  dominant  factor. 

The  slope  of  the  Rayleigh  scattering  curve  in  Figure  40  be¬ 
tween  27000  and  28000  cm'^  is  within  10%  of  the  slope  in  the 
efficiency  curve,  Figure  37,  for  the  same  region. 

Beam  Radius.  The  total  efficiency  of  the  laser  system 
scales  as  one  over  the  square  of  the  laser  beam  radius  at 
the  target,  therefore,  beam  radius  determination  should  be 
very  significant.  As  in  the  case  of  atmospheric  transmission, 
values  of  the  individual  components  of  the  beam  radius  were 
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output  from  the  main  program  for  examination. 

The  beam  radius  which  results  from  beam  diffraction 
effects  was  expected  to  scale  inversely  with  the  wavenumber 
of  the  laser  radiation,  favoring  higher  wavenumber.  Exami¬ 
nation  of  the  plot  in  Figure  41  indicates  that  this  was 
indeed  the  case  and  hand  calculations  confirmed  the  values. 
Further,  the  effects  of  focusing  can  be  seen  as  the  beam 
radius  at  the  target  is,  for  high  wavenumbers,  reduced  to 
values  smaller  than  the  laser  output  optics.  Although  not 
plotted,  the  effects  of  jitter  produced  the  expected  constant 
radius  across  the  spectrum. 

For  this  analysis,  we  purposely  selected  atmospheric 
conditions  which  were  unfavorable  for  turbulence.  This 
choice  proved  to  be  very  significant.  The  beam  radius  which 
results  from  turbulence  effects  was  expected  to  be  a  very 
weak  function  of  wavenumber  (to  the  1/5  power)  and  to  favor 
the  lower  wavenumbers.  Figure  42  shows  this  weak  wavenumber 
dependence.  One  might  reasonably  expect  that  the  stronger 
wavenumber  dependence  in  the  diffraction  effect  to  win  out 
over  the  weaker  dependence  in  turbulence.  However,  this 
proved  not  to  be  the  case.  Even  with  our  85%  reduction  in 
turbulence  induced  beam  radius  to  account  for  some  degree  of 
atmospheric  compensation,  this  radius  was,  in  the  low  wave- 
number  end  of  the  spectrum,  twice  the  diffraction  radius  and 
fully  10  times  the  diffraction  radius  in  the  high  wavenumber 
end  of  the  spectrum.  The  result  being  that,  when  the  indi¬ 
vidual  beam  effects  are  combined,  the  total  beam  radius, 


rather  than  decreasing,  increases  slightly  with  rising  wave- 
numbers.  Likewise,  the  overall  system  efficiency  decreases. 
The  total  beam  radius  is  depicted  in  Figure  43.  It  was  clear 
that  with  the  original  choice  of  parameters,  turbulence 
induced  effects  dominated  the  total  beam  radius. 

FEL  efficiency  and  target  efficiency  were  not  examined 
separately  since  their  values  did  not  depend  upon  the  selec¬ 
tion  of  baseline  parameters.  From  Figure  11,  it  can  be  seen 
that  the  FEL  efficiency  is  a  weak  function  of  wavenumber, 
favoring  lower  wavenumbers.  On  the  other  hand,  the  target 
efficiency  scales  somewhat  more  strongly  and  favors  the 
higher  wavenumbers.  Target  efficiency  ranges  from  a  value 
of  1.0  to  .55  over  the  selected  wavenumber  region. 

When  just  the  effects  of  aerosols,  ozone  absorption, 
Rayleigh  scattering,  beam  radius,  FEL  efficiency,  and  target 
efficiency  are  combined,  a  curve  is  formed  which  perfectly 
duplicates  the  top  line  structure  of  the  efficiency  curves 
computed  by  the  model.  This  curve  is  plotted  in  Figure  44. 

In  its  compressed  wavenumber  scale,  the  optimum  peak  is 
clearly  more  discernable,  although  no  more  valid.  Also  more 
discernable  is  the  weak  ozone  attenuation  which  proves  to  be 
a  dominating  factor  in  the  15000  to  18200  era"'*'  region. 
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Extensions  to  the  Base  Case 

With  the  representative  parameters  chosen  for  the  laser 
system,  we  were  able  to  define  an  optimum  wavenumber  value. 
However,  within  a  relatively  broad  wavenumber  region  either 
side  of  this  optimum  value  the  system  efficiency  was  approx¬ 
imately  equal  and  the  wavelength  selection  essentially  in¬ 
different.  Furthermore,  the  assumptions  and  inaccuracies 
inherent  to  the  model  tend  to  make  the  concept  of  an  opti¬ 
mum  wavenumber  even  more  ambiguous.  Nonetheless,  the  con¬ 
cept  is  useful  in  examining  the  sensitivities  in  the  wave¬ 
length  selection  process  for  the  free  electron  laser.  In 
the  remainder  of  this  report,  we  will  concentrate  on  how  the 
optimum  wavenumber  departs  from  the  value  found  in  the  base¬ 
line  analysis  when  certain  controllable  and  uncontrollable 
parameters  are  altered. 

To  do  this,  the  analysis  process  will  be  greatly  simpli 
fied.  First,  thermal  blooming  effects  will  be  ignored. 
Further,  the  molecular  absorption  data  from  FASCODE  will  be 
excluded.  By  doing  so,  the  overall  system  efficiency  become 
a  smooth  function  of  wavenumber,  resembling  the  curve  in 
Figure  44.  Through  such  an  analysis,  the  approximate  loca¬ 
tion  of  the  peak  can  be  determined.  Undoubtedly,  this  point 
may,  at  times,  fall  on  a  molecular  absorbing  feature;  yet, 
the  purpose  is  not  to  determine  an  exact  wavenumber  for  FEL 
operation  but  rather  to  examine  the  degree  of  impact  that 
certain  parameters  have  on  fixing  the  spectral  region  where 
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such  a  wavenumber  can  be  found. 

Using  this  simplified  approach,  the  baseline  case  was 
recalculated.  This  time,  in  the  absence  of  molecular  absorp 
tion  data,  the  peak  wavenumber  occurred  at  14405.0  cm*'*'  with 
a  computed  efficiency  of  .5085%.  This  wavenumber  does  in¬ 
deed  fall  on  an  absorption  feature,  and  explains  the  slight 
shift  in  the  peak  wavenumber  (~  58  cm  ^ )  and  efficiency  from 
the  results  from  the  full  model. 

Effect  of  Turbulence  Compensation.  As  discussed  above, 
the  effects  of  turbulence  actually  caused  the  beam  radius 
to  increase  with  increasing  wavenumbers  in  the  baseline 
case,  rather  than  decrease.  Therefore,  it  is  expected  that 
turbulence  effects  would  have  a  profound  influence  upon  the 
peak  wavenumber  along  the  efficiency  curve.  An  experiment 
was  conducted  with  the  simplified  model,  maintaining  all 
laser  parameters  equal  to  the  baseline  case  except  for  the 
scaling  term  for  atmospheric  compensation.  This  was  varied 
from  .75  to  1.00  in  intervals  of  .025.  The  results  of  this 
analysis  can  be  found  in  Table  I. 

It  can  be  seen  that  altering  the  degree  of  atmospheric 
compensation,  drastically  affects  the  peak  in  the  efficiency 
curve;  nearly  the  entire  spectrum  under  study  is  covered  by 
this  experiment.  However,  even  more  significant  are  the 
order  of  magnitude  increases  in  overall  system  efficiency 
one  achieves  with  enhanced  compensation.  When  the  beam 
radius  approaches  its  diffraction  limited  value,  both 


TABLE  I  -  Variation  of  Turbulence  Compensation 


Peak 

Peak 

Wavenumber 

Efficiency 

w  cm"  ^ 

n  7. 

9435 

.1333 

12910 

.1728 

14005 

.  2335 

14390 

.  3324 

14405 

.  5085 

14405 

.  5836 

20200 

1 .8306 

22200 

3.9726 

23010 

5 .8346 

27870 

24 . 7141 

Degree  of 
Compensa  tion 


Mi] 


Wavenumber 

Shift  from 
Baseline  Case 


-4970 

-1495 

-400 

-15 

0 

0 

+  5  7  9  5 
+  7  795 


+  8605 


65 
+14590 


TABLE  II  -  Variation  of  Turbulence  Compensation 
(No  Ozone  Absorption) 


Degree  of 
Compensation 


Peak 

Wavenumber 

X)  -1 


9435 
12910 
14005 
14390 
14405 
1  52  70 
16800 
18185 

23010 
27  8  7  0 


WHIUi] 


Peak 

Efficiency 
b  7. 


.1333 
.1728 
-.2335 
.3412 
.5125 
.  9233 
1.3812 
2.8573 

5 .8346 
24.7141 

-  1  : 


Wavenumber 
Shift  from 
Baseline  Case 


497 

149 

400 

15 

0 

865 

395 


+  3  780 
+  8605 


efficiency  and  peak  wavenumber  increase  dramatically.  A 
plot  of  the  results  for  100%  turbulence  compensation  is 
found  in  Figure  45  to  illustrate  this  point. 

Admittedly,  these  results  are  somewhat  artificial  for 
two  reasons:  1)  the  degree  of  compensation  does  not  take 
into  account  the  fact  that  limits  do  exist  as  to  how  much 
correction  can  actually  be  accomplished  through  adaptive 
optics  on  a  real  laser  beam,  and  2)  we  have  left  the  con¬ 
tribution  from  ozone  absorption  in  the  calculations,  causi. 
the  efficiency  peaks  to  remain  on  either  side  of  the  weak 
ozone  absorption  feature  in  the  mid-visible.  To  address 
the  latter  inconsistency,  a  second  experiment  was  accom¬ 
plished.  This  time,  the  contributions  from  ozone  attenua¬ 
tion  were  excluded.  The  results  of  that  analysis  can  be 
found  in  Table  II. 

Without  the  ozone  contributions,  the  baseline  peak 
wavenumber  remains  the  same,  yet,  the  system  efficiency 
changes  slightly.  In  this  experiement,  peak  values  of  sys- 
tme  efficiency  shifted  at  somewhat  more  regular  intervals, 
yet  the  results  and  conclusions  remain  essentially  the  same 

Effects  of  Altitude  and  Aerosols.  The  height  above  se 
level  of  the  laser  and  the  selected  aerosol  model  were  ana¬ 
lyzed  together.  The  baseline  case  assumes  a  'best  case' 
aerosol  model  based  on  a  rural  environment  and  high  visibil 
ity  (23  km)  conditions.  The  height  of  the  laser  in  the  bas 
line  case  (1.0  km)  was  admittedly  pessimistic.  In  this 
section  we  will  examine  the  effects  of  increasing  the  laser 
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altitude  and  decreasing  the  visibility  due  to  aerosol  con¬ 
stituents  in  the  atmosphere. 

Sea  level  visibilities  of  5,  15,  and  23  km  were  chosen 
for  the  aerosol  model  in  LOWTRAN  while  the  Rural  aspect  of 
the  model  in  the  boundary  layer  was  retained.  The  height  of 
the  laser  was  increased  in  1.0  km  increments  to  3.0  km.  The 
results  of  this  third  experiment  are  found  in  Table  III. 

As  expected,  both  the  peak  wavenumber  and  the  peak 
efficiency  decrease  as  the  visibility  decreases  when  the 
laser  remained  at  1.0  km.  The  decrease  is  slight  between 
23  km  and  15  km  visibilities,  but  becomes  significant  be¬ 
tween  15  km  and  5  km. 

At  2.0  km  and  above,  the  laser  is  out  of  LOWTRAN 's 
boundary  layer  model,  and  hence,  the  differences  in  sea 
level  visibility  have  no  effect.  This  is,  of  course,  not 
realistic.  Even  at  3.0  km.  surface  effects  will  determine 
actual  aerosol  concentrations.  Unfortunately,  LOWTRAN 
assumes  that  the  surface  is  at  sea  level  even  though  the 
initial  altitude  for  its  calculations  might  be  3.0  km.  The 
user  can  work  around  this  problem  by  inserting  his  own 
aerosol  model.  However,  these  results  can  illustrate  cer¬ 
tain  points.  Raising  the  height  of  the  laser  from  1.0  to 
2.0  km  results  in  a  dramatic  shift  in  the  peak  wavenumber 
to  higher  wavenumbers.  Increasing  the  altitude  by  another 
kilometer,  shifts  that  peak  still  farther,  even  though  the 
aerosol  attenuation  is  approximately  the  same.  This  is,  of 


TABLE  III  -  Variation  of  Altitude  and  Sea  Level  Visibility 


Altitude 


5.0  km 


Sea  Level  Visibility 

15.0km  23.0km 


Peak  Pea 


1.0  km  12930  .4454  14305 
2.0  km  22200  .6535  22200 
3.0  km  23910  .7147  23910 


v  cm  r)  7. 


14405  .5085 

22200  .6535 

23910  .7147 


TABLE  IV  -  Variation  of  Turbulence  Compensation 
(  Laser  at  3.0  km  Altitude  ) 


Degree  of 
Compens a  t ion 


Peak 

Wavenumber 

„  -1 
v  cm 

Peak 

Efficiency 

n  7. 

Wavenumber 

Shift  From 

Baseline  Case 

2  3645 

.1811 

-  265 

23850 

.  3208 

-60 

23910 

.7147 

0 

24505 

1.4328 

+  595 

2  5600 

4.1161 

+  1690 

29000 

+  5090 

course,  due  to  the  much  decreased  Rayleigh  scattering  con¬ 
tribution.  The  atmospheric  density  is  down  by  a  factor  of 
3/4  from  the  sea  level  value  at  an  altitude  of  just  3.0  km 
and  the  strong  wavenumber  dependence  in  Rayleigh  scattering 
causes  a  significant  shift  in  the  efficiency  peak.  The  com¬ 
puted  efficiency  curve  for  a  3.0  km  laser  altitude  is 
depicted  in  Figure  46. 

Effects  of  Turbulence  Compensation  -  (2).  The  signifi¬ 
cant  impact  of  the  height  of  the  laser  on  the  peak  wavenumber 
suggested  a  re-examination  of  the  effects  of  turbulence  com¬ 
pensation  when  the  laser  is  raised  to  a  much  higher  altitude. 
An  experiment  was  accomplished  with  the  laser  at  an  altitude 
of  3.0  km  and  with  a  23  km  visibility,  Rural  aerosol  model. 
The  turbulence  compensation  factor  was  varied  from  .75  to 
1.00  in  intervals  of  .05.  The  results  of  this  experiment 
are  found  in  Table  IV. 

The  effect  of  turbulence  compensation  on  the  efficiency 
of  the  system  is  still  great,  changing  two  orders  of  magni¬ 
tude  over  the  selected  range  of  the  compensation  parameter. 
Although,  compensation  is  still  quite  important  in  deter¬ 
mining  the  peak  wavenumber,  the  impact  is  much  less  at  the 
higher  laser  altitude.  Here,  enhanced  aerosol  and  Rayleigh 
scattering  transmission  tend  to  play  a  much  greater  role  in 
determining  the  overall  weapon  system  efficiency. 


Effects  of  Aperture  Size.  The  radius  of  the  laser  out¬ 
put  optics  was  probably  under  estimated  in  the  baseline  case; 


.4890 
.  9492 
1.5121 


5 . 5 


14305 


2 . 1444 


-100 


With  larger  optics,  the  same  diffraction  angle  of  the 
laser  beam  can  be  achieved  through  longer  wavelengths 
(lower  wavenumbers);  therefore,  as  the  results  show,  increas¬ 
ing  the  size  of  the  optics  decreases  slightly  the  peak  wave- 
number.  However,  the  effect  on  the  value  of  the  peak 
wavenumber  is  very  small.  A  more  significant  impact  is 
found  on  the  overall  efficiency  of  the  weapon  system. 


VII.  Conclusions  and  Recommendations 


Unfortunately,  time  precluded  a  more  thorough  examina¬ 
tion  of  the  factors  which  impact  the  selection  of  an  optimum 
wavelength.  This  chapter,  however,  will  summarize  those  con 
elusions  we  have  made  during  the  study  and  briefly  offer 
recommendations  for  further  research. 

Conclusions 


Many  wavelength  dependent  relationships  which  might 
impact  upon  the  selection  of  an  operating  wavelength  for  a 
ground  based  free  electron  laser  weapon  system,  were  consid¬ 
ered  in  this  study.  Thermal  blooming  was  found  not  to  be  a 
factor  in  the  selection  process.  During  the  analysis,  no 
single  wavelength,  nor  even  a  small  range  of  wavelengths, 
could  be  unambiguously  identified  as  optimum  for  the  strate¬ 
gic  mission.  For  the  system  parameters  we  chose,  an  optimum 
wavelength  was  determined  to  be  0.697  micrometers  (14347 
cm  ^ ) .  However,  this  result  is  inconclusive  for  a  number 
of  reasons. 

System  efficiencies  within  only  1%  of  the  optimum 
value  could  be  found  over  a  3000-4000  cm"*'  wavenumber  re¬ 
gion  surrounding  the  peak  wavenumber.  Taking  into  account 
the  errors  and  inaccuracies  in  the  model,  which  we  have 
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estimated  to  be  5%,  the  optimum  wavenumber  could  fall  any¬ 
where  between  11000  cm  ^  (0.91  micrometers)  and  24000  cm"''' 
(0.42  micrometers). 

In  addition,  it  is  expected  that  the  laser  system 
parameters  which  were  chosen  for  the  analysis,  could  be 
unrealistic  or  inappropriate  for  a  specific  mission.  Real¬ 
izing  this  fact,  we  sought  to  determine  how  sensitive  our 
optimum  wavenumber  value  was  to  changing  parameters. 

When  the  degree  of  compensation  for  atmospheric  tur¬ 
bulence  is  varied,  the  peak  in  the  efficiency  curve  shifts 
significantly.  Even  with  the  excellent  compensation  factor 
we  chose  for  the  analysis,  the  turbulence  induced  beam 
radius  dominates  the  total  beam  radius  and  a  modest  change 
in  the  ability  of  an  adaptive  optics  system  to  correct  for 
turbulence  can  have  a  dramatic  influence  on  the  selection  of 
a  wavelength.  Turbulence  has  the  most  profound  effect  on 
the  overall  system  efficiency  which  can  change  two  orders 
of  magnitude  depending  on  the  degree  of  compensation.  The 
major  conclusion  which  can  be  drawn  from  this  discussion  is 
already  well  known.  Unless  one  does  an  excellent  job  of 
correcting  for  atmospheric  turbulence,  the  beam  radius,  and 
hence  the  system  efficiency,  becomes  quickly  dominated  by 
turbulence  effects  and  little  advantage  accrues  from  having 
shorter  wavelengths,  larger  apertures,  focusing  and  high 
beam  quality.  We  would  add  that  the  degree  of  compensation 
has  a  direct  and  significant  bearing  on  the  selection  of 
laser  wavelength. 


When  the  mean  sea  leveL  altitude  of  the  laser  is  raised, 
the  peak  wavenumber  increases  dramatically  due  to  the  greatly 
reduced  aerosol  concentrations  and  Rayleigh  scattering  ef¬ 
fects.  At  an  altitude  of  3.0  km,  the  optimum  wavelength 
shifts  to  approximately  0.41  micrometers.  Raising  the  al¬ 
titude  of  the  laser  produces  a  modest  increase  in  efficiency 
as  well.  At  the  higher  altitudes,,  the  impact  of  turbulence 
compensation  on  the  peak  wavenumber  lessened  but  still 
significant.  However,  the  ability  to  correct  for  turbulence 
remains  just  as  important  at  the  higher  altitudes  to  the 
overall  weapon  system  efficiency. 

Increasing  the  size  of  the  output  optics  will  greatly 
enhance  the  system  efficiency,  yet  has  little  impact  on 
the  optimum  wavenumber  for  laser  operation. 

In  summary,  laser  altitude  is  a  strong  influence  on 
the  peak  wavelength  while  turbulence  compensation  also  pro¬ 
duces  large  shifts  in  peak  wavelengths  but  enormous  changes 
in  the  overall  system  efficiency. 

Recommendations 

This  analysis  has  hopefully  examined  most  of  the  factors 
bearing  on  the  problem;  however,  the  actual  selection  pro¬ 
cess  will  likely  be  determined  by  further  study. 

The  amount  of  correction  for  atmospheric  turbulence 
which  is  eventually  demonstrated  will  impact  heavily  on  the 
selection  decision.  In  fact,  the  peak  wavelength  may  be 


totally  determined  by  the  wavelength  dependence  in  an 
adaptive  optics  system.  In  retrospect,  there  should  have 
been  more  of  an  attempt  to  include  realistic  adaptive  optics 
capabilities  in  this  analysis. 

In  addition,  the  efficiencies  of  optical  devices,  such 
as  cavity  mirrors  and  focusing  elements,  ignored  by  this 
study,  will  likely  have  an  important  influence  once  a  nomi¬ 
nal  wavelength  is  found  and  the  optical  materials  selected. 

Lastly,  the  wavelength  dependence  in  the  target  coup¬ 
ling  efficiency  remains  a  very  open  question.  It  will 
probably  be  eventually  answered  through  experimental  studies 
rather  than  theoretically  scaling  relationships.  The  pro¬ 
gram  for  this  analysis  has  been  constructed  in  such  a  way 
as  to  easily  allow  the  addition  of  specific  spectral  target 
data  at  a  later  time. 
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Modifications  to  the  Code  -  LOWTRAN 


This  appendix  describes  the  modifications  which  were 
accomplished  during  this  study  on  the  code  LOWTRAN,  version 
6.  The  purpose  was  two-fold:  1)  to  extract  absorption  and 
scattering  optical  depths  for  each  of  the  32  altitude  lay¬ 
ers  in  LOWTRAN,  and  2)  to  extract  information  on  the 
atmospheric  model  for  use  in  the  main  program.  These  mod¬ 
ifications  are  listed  as  follows: 

A.  In  the  main  program,  LWTRN6,  two  files,  TAPE8  and 
TAPE15,  are  declared  in  the  PROGRAM  statement,  line  110. 

B.  In  subroutine  STDMDL,  input/output  unit  15  is 
opened  between  lines  250  and  255,  and  designated  the  file 
name,  MODEL.  Into  this  file,  the  atmospheric  data  will  be 
written.  Later  in  this  subroutine,  between  lines  650  and 
655,  a  WRITE  statement  is  inserted  to  extract  the  atmospheric 
pressure,  temperature,  and  altitude  of  the  33  altitude 
boundaries  in  LOWTRAN.  An  implied  DO  Loop  structured  is 
utilized . 


WRITE  (15)  (  PM<  I) ,  TM(I) ,  ZM(I),  1=1,33  ) 


C.  After  LOWTRAN  computes  the  path  geometry  through 
the  atmosphere,  the  average  values  of  pressure,  temperature 
and  density  for  each  of  the  32  altitude  layers  are  also 
written  to  the  file  MODEL.  This  is  accomplished  in  sub¬ 
routine  RFPATH.  This  subroutine  loops  through  each  layer 
and  computes  these  values  but  does  not  store  them.  There¬ 
fore,  it  is  necessary  to  create  arrays  in  which  to  store 
these  values.  This  is  done  in  the  DIMENSION  statement  in 
line  200  with  the  addition  of  variables  PBAR  (35),  TBAR(35) 
and  RHOBAR(35).  Within  the  layer  loop,  arrayed  variables 
are  substituted  for  the  single  dimension  variables  of 
PBAR,  TBAR,  and  RHOBAR.  This  is  done  at  lines  560-570  and 
at  lines  860-870.  Finally,  outside  of  the  layer  loop, 

the  arrays  are  written  to  the  file  MODEL  after  line  915. 
The  file  is  closed  immediately  following  this  operation. 

WRITE  (15)  (  PBAR(I),  TBAR(I),  RHQBAR(I) ,  1=1.32  ) 

D.  The  most  convenient  place  to  extract  absorption 
and  scattering  data  is  in  subroutine  TRANS.  This  is  accom¬ 
plished  by  first  accessing  a  pre-existing  common  block  of 
data  through  a  COMMON  statement. 

COmON  /  RFRPTH  /  ZP(35)  ,PP(35)  .TP(35)  ,RFNDXP<35)  ,3P(35) . 

1  PPSHK35)  ,TPS(JM(35)  ,RH0SUM(35)  ,DENP(  16.35)  ,AMTP(  16.35) 


The  arrayed  variable  AMTP(X,Y)  from  above  represents  the 
column  density  of  the  atmospheric  constituent  X  over  the 
altitude  layer  Y.  Next,  the  arrayed  variables  RSCAT(35), 
AEXT  (35),  ABBS ( 35 ) ,  CONT(35),  and  OZONE(35)  are  dimen¬ 
sioned,  and  the  input/output  unit  8  is  opened  and  given 
the  file  name,  LTDATA.  These  variables  represent  the  con¬ 
tributions  of  Rayleigh  scattering,  aerosol  scattering, 
aerosal  absorption,  continuum  absorption,  and  ozone  absorp¬ 
tion.  The  arrayed  variable  ABB(X)  from  LOWTRAN  corresponds 
to  the  cross-sectional  probability  for  constituent  X.  The 
following  is  then  inserted  between  lines  1250  and  1255. 

DO  359  M»l,32 

RSCAT(N)=  AB8(8)  *  AMTP(12,N) 

A=  XXI  if  AMTP(7,N)  +  XX2  i  MTP(12.N) 

B=  XX3  *  AMTP(  13, N)  +  XX4  if  f¥fTP(J4,N) 

AEXT(N)=  A  +  B 

A=  YY1  if  AMTP(7,N)  +  YY2  if  f¥fTP(12.N) 

8=  YY3  i  MTP<I3,N)  ♦  YY4  if  fWP(14,N) 

ABBS(N)=  A  +  B 

C0NT(N)=  ABB (5)  if  tffTP<5,N)  ♦  ABB(?)  if  tffTPfP.N)  + 

1  ABB(  18)  i  AMTP( 18  ,N) 

OZONE (N)=  ABB (3)  if  tf1TP(8,N) 

358  CONTINUE 

WRITE  (8)  0,  (RSCAT(l) , AEXT  Cl) ,AB8S(I) ,020NEC I) . 

1  CONT(I),  1=1,32) 

Outside  of  the  wavenumber  loop  in  subroutine  TRANS,  after 
line  1565,  the  unit  8  can  be  closed. 
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Modifications  to  the  Code  -  FASCODE 


This  appendix  describes  the  modifications  on  the  code 
FASCODE  which  were  accomplished  to  extract  molecular  absorp¬ 
tion  data  for  each  altitude  layer.  These  modifications  are 
as  follows: 

A.  In  the  main  program,  FASC0D1,  an  additional  file, 
TAPE14,  is  declared  in  the  PROGRAM  statement,  line  110. 

Also,  this  input/output  file  is  opened  and  assigned  the  name 
FCDATA,  following  line  1362  in  the  main  program.  An  ENDFILE 
statement  and  a  CLOSE ' s tatement  on  this  file  are  then  in¬ 
serted  following  line  1470. 

B.  Subroutine  LAYER  controls  the  program  looping  over 
the  altitude  layers.  In  order  to  distinguish  these  differ¬ 
ent  layers  in  the  data  on  file  FCDATA,  it  was  necessary  to 
place  a  marker  after  each  layer  calculation.  This  is  accom¬ 
plished  by  inserting  between  lines  5020  and  5025  the  follow¬ 
ing. 


LL  =  0 
XX  =  0.0 

WRITE  (14)  XX, XX, 11 


C.  In  subroutine  PATH,  the  sampling  interval  for 
FASCODE  is  computed.  We  wished  to  fix  this  interval  to  0.2 


wavenumbers.  Therefore,  line  7530  of  this  subroutine  is 
replaced  by: 


38  OU  =  0.2 

D.  In  subroutine  PANEL,  line  16470  is  modified  to  pass 
the  number  of  the  current  layer  to  subroutine  R1PRNT.  One 
must  remember  to  specify  a  value  greater  than  zero  for  the 
parameter  NPTS  on  the  FASCODE  input  cards  in  order  for  this 
subroutine  to  be  called. 


IF  (NPTS  .GT.  9)  CALL  R1PRINT  (R1.NL0.KFILE, LAYER) 

E.  Finally,  subroutine  R1PRNT  is  re-written  entirely  and 
is  presented  here  in  its  entirety.  The  variable  R1  is  the  op¬ 
tical  depth  for  wavenumber  VJ  over  altitude  layer  NUM. 


SUBROUTINE  R1PRNT  (R1.JL0.MFILE.NUM) 

C(W«3N/XPfmA>  IP  .U2P ,  OMP  ,NLIM  .NSHFT.NPTS 
C0t1CN/IFIl/IRD,IPR,IPU.NQPR 

DIMENSION  Rl(l) 

JHI=JL0  ♦  NLIM  -  1 
00  18  KK=1,NLIM 
J=JL0  +  KK  -1 

UJ=U1P  +  FLOAT (J-JLO)  X  DYP 
WRITE  (14)  UJ,  Rl(J),  NflJM 
18  CONTINUE 
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Proaram  Listina  -  S 


This  program  is  designed  to  sort  the  molecular  absorp¬ 
tion  data  generated  by  FASCODE  into  a  format  which  is  compat 
ible  to  the  main  program  and  LOWTRAN.  The  program  reads  a 
block  of  2000  wavenumbers  for  a  given  altitude  layer,  and 
then  skips  via  the  subroutine  POST  to  the  appropriate  point 
in  the  next  altitude  layer.  When  all  32  layers  have  been 
read,  the  data  is  written  to  file  and  the  process  repeats 
with  the  next  2000  wavenumbers.  The  output  on  file  SDATA 
now  contains  molecular  absorption  data  over  all  layers  for 
a  single  wavenumber.  The  program  is  as  follows: 


PROGRAM  SORT 

REAL  M0LABS ( 2080 , 32)  ,WAWUM(20Ue>  ,V,AB8 
INTEGER  1. LEVEL. KCOUNT 

0PEN( 14,FILE='FCDATA') 

OPENt 12,FIL£=/S0ATA/) 

KCOJNT=0 


10  1=1+1 

DO  100  K= 1,2000 

READ  (14,END=500)  V, ABB, LEVEL 


IF  (LEVEL  .EQ.  3)  THEN 
IF  (I  .LT.  32)  THEN 
IF  (KCOUNT  .EQ.  i)  THEN 
GO  TO  10 
ELSE 

CALL  POST  (KCOUNT) 

GO  TO  18 
END  IF 
ELSE 

STOP  1889 
END  IF 
ELSE 

IF  (I  .NE.  LEVEL)  THEN 
STOP  2988 
END  IF 
HAWUM(K)=V 
MOLA8S(K, LEVEL) =A8B 
END  IF 
180  CONTINUE 

200  READ( 14,END=250)  V, ABB, LEVEL 
IF  (LEVEL  .NE.  8)  TfCN 
GO  TO  200 
ELSE 

IF  (I  .LT.  32)  THEN 
IF  (KCOUNT  .EQ.  9)  THEN 
GO  TO  10 

asE 

CALL  POST  (KCOUNT) 

GO  TO  10 
END  IF 

asE 

GO  TO  258 
END  IF 
END  IF 

258  DO  388  N=1,K-1 

WRITE  (12)  WAWNUM(N) ,(MOLABS(N,I) ,  1=1,32) 
309  CONTINUE 

KCOUNMCOUNT+1 
REMIND  14 
CALL  POST (KCOUNT) 

1=8 

GO  TO  10 

588  DO  498  N»1,K-1 

WRITE  (12)  l<ttWUM<N),  (MOLABS(N.I) ,  1=1.32) 
688  CONTINUE 
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APPENDIX  D 


Program  Listing  -  WAVLEN 


PROGRAM  HAWLEN  ( INPUT, OUTPUT, TAPE5=  INPUT ,TAPE6=0UTPUT) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  PROGRAM  BY 

X  CAPT  DOUGLAS  E.  KGHLHEPP 
X  AFIT/ENS  ,  GS0-83D 

X  WRIGHT-PATTERSON  AFB,  OHIO 

X 

X  THIS  PROGRAM  IS  DESIGNED  TO  FIND  THE  OPTIMUM  WAVELENGTH 
X  FOR  A  GROUND-BASED  FREE  ELECTRCN  LASER 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  M0LABS(32) ,AEXT(32) ,ABBS<32) .RSCATC32) ,OZONE<32) . 

1  ATTEN(33)  ,A8S0RB(33)  ,DELZ<32)  ,P(33)  ,T(33)  ,ALT(33) . 

2  PBARI32) ,T8AR(32) ,RHOBAR(32) ,IAVG,PBEST,VBEST, JITTER, HOL, 

3  F0CRNG,LAM3A,HABG,If1SL,NI,IREL.ILASER,C0NT(32)  .RANGE, 

4  COMPEN , EBEST , VPBEST , EFF , ETA , PROPEFF 
INTEGER  H1.IFLAG.NSTEP 

COtlON/  OPOPTH/M&ABS ,  AEXT  .ABBS , RSCAT , OZONE .  CONT 

COimi/ATM/P,T,ALT,PBAR,TBAR,RHOBAR 

COWON/RAD/HOL  .JITTER ,  BEANO ,  FOCRNG ,  RHOTURB  ,RB ,  COMPEN 

DATA  (DELZ(I) .1=1 .253/25X 1900 .9/ 

DATA  (DEL2(I) ,I=26,30)/5X5880.8/ 

DEL2(  31)  =20000.0 
DEL2( 32) =30000.0 

X  CONSTANTS  AND  STARTING  VALUES 
PI=3.141<4 
IFLAG=0 
EBEST=0.0 
:3EST=0.0 
ATTEN(33)=0.0 
ABSORB (33) =0.0 

X  OPEN  DATA  FILES 

0PEN(8,FILE='IT0ATA') 

OPEN( 14,FILE='FCDATA') 

OPENC 15,FILE="M0DEL') 

OPEN( 10,FILE='EFFIC') 

REWIND  3 


161 


•j  vvvvSV'.V' 


w, 


•» 


* 


REWIND  14 
REWIND  15 

XREAD  IN  PARAMETERS  FOR  ANALYSIS 
XREAD  CARD  #1 

READ(5,  108) AVGPWR.R0, VI, V2, RANGE 
108  FORMAT  (5F12.2) 

RANGE=RANGEX 1800.0 

XfiEAO  CARO  #2 

REAO(5.101)  FOCRNG.BEAMQ. JITTER. HOL.COMPEN 
101  FORMAT  (5F12.2) 

FQCRNG=FOCRNGX 1000.0 
Hl=NINT(HOL)  +  1 
HOL=HOU  1000.8 

IF  (yi.GT.02)  THEN 

PRINT' ('WAVENUMBER  PARAMETER  ERROR1) ' 

STOP  1008 
END  IF 

it  READ  IN  ATMOSPHERIC  DATA 
CALL  MDLATN 

X  BEGIN  COMPUTATION  OF  YURA'S  COHERENCE  DIAMETER  FOR  TURBULENCE 
X  SELECT  AN  AVERAGE  WAVELENGTH  OF  .6  MICRONS 
UMD*=.<iE-4 

XCQNSTANTS  AND  INITIAL  VALUES 
HABG*0.0 
SUM=0.0 
DELH=50.0 
0=247.4 
USTAR=0.35 

X  BEGIN  SIM1ATI0N  NEAR  EARTH'S  SURFACE  WITH  WYNGMRO'S  MODEL 

DO  110  L=50, 950.58 

HABG=L+(DELH/2.0) 

UTSL=HOL  +  HA8G 

TEMP*  T(H1)  +  (  (HABG/DELZCHl))  X  (T(H1+1)  -  T(HD)  ) 

SCALE*  P(H1+1)  X  T(Hl)/(  P(H1)  X  T(H1+1)  ) 

PR£S=  TEMP  X  (  P(H1)/T(HD)  X  SCALE  XX  (  HABG/DEL2(H1)  ) 

DENOM*  (14000.0  X(USTAR  XX  3.0))/<HABGXQ) 

OENO*  (DENOM  ♦  1.8)  XX  .667 


CT2=  (2.0E-3  X  (Q  X*  1.333)  *  (WiBG  XX  (-1.333) ,>)/DENCM 
0*2-  CT2  X  (DN0T(FRE3,TEMP.L^IDh)  XX  2.3) 


A=(  (RANGE-HABG) /RANGE)  XX  1.667 
A=  A  X  CN2  X  OELX 

SUH-SUM  +  A 

118  CONTINUE 

DO  115  MM*  1,25 

DO  113  N— 1 ,2 

ZBAR=ALT(I)  +  (2.8XN  -1.0)/<DELZ(I)X4.0)  -  HOL 

CN2=  2.23K1.8E-5  X  18.0XX(-.3)  X  ZBAR)  XX10.0  X  EXP(-ZBAfi/1080.8) 
CN2=  CM2  ♦  1.0E-14  X  EXP(-ZBAR/ 1500.0) 

CN2=  CN2  X  2.7 

ApuRANGE-ZBAR) /RANGE)  XX  1.447 
Af*  CN2  X  A  X  DELZ(I) 

SUM*SUM  ♦  A 

113  CONTINUE 

115  CONTINUE 

WRITE  (4,117)  SUM 
117  FORMAT  ('  SUM*  '.E12.4) 

X  BEGIN  WAUENWBER  LOOP 

X  READ  IN  FASCQDE  DATA 

120  IF  (IFLAG  .EQ.  0)  THEN 

CALL  FCREAO  (0.  IFLAG. 01, 02) 

END  IF 

IF  (IFLAG  .EQ.  1)  THEN 
0=0+ 1.0 
DO  125  1=1,32 
MOLABS(I)=8.0 
125  CONTINUE 
END  IF 

IF  (V  .GT.  02)  THEN 
GO  TO  388 
END  IF 


XREAO  IN  LOWTRAN  DATA 


CALL  lTREAO  (V,  I  FLAG) 

IF  4. 1  FLAG  .EG.  2)  THEN 
GO  TO  380 
END  IF 

i  COWERT  WAVENU46ER  TO  WAVELENGTH  IN  METERS 
LAMDAf1.8E-2/V 

i  COMPUTE  TOTAL  TRANSMISSION  ANO  THERMAL  BLOOMING  IN  THE  FIRST 
*  LAYER  TO  DETERMINE  IF  IT  IS  WORTHWHILE  PROCEEDING 
TRAN=8.8 
DO  158  I*H1.32 

ATTENC  I)=M0LA8S(  I)  +AEXT(  I)  tAB8S(  l)  +RSCAT(  I)  +Q20NEC  I)  +CONT  ( I) 
ABSORB!  I)=MOLABS(  I)  +AB8S(  I)  +020NEC I) 

TRAN=TRM+ATTEN(I) 


158  CONTINUE 

TftW=EXP(-TRM) 

NI=-0NDT(PBAR(H1) .TBAR(Hl) ,LAMDA)  *  AVGPWR 
NI=NI*A8S0RB(H1)*DELZ<H1)*DELZ(H1) 

NI=NI/(PIIRH0BAR(H1)  *988 .8*VX(H0L)  *R8*R0*R8) 

IF  (TRAN  .LE.  8.18  .OR.  NI  .GE.  38.8)  THEN 
EFF=8.8 

WRITEU8)  V.EFF 
GO  TO  128 
ENO  IF 

IFINISH  WAVELENGTH  OEPENOENT  CALCULATIONS  OF  YURA'S  CCKERFNCE  DIAMETER 
RHOTURB*<  5 . 88*P  I  *  PI  *  SUM)  /  ( LAMOAifLAMDA) 

RHOTURB=RHOTURB  it  (-.68) 

4 START  THERMAL  BLOOMING  CALCULATIONS 

NSTEP=4 

NI=8.8 

4  START  OUTER  LOOP 
DO  358  I*Hl,32 
0*8.8 

22=ALT(I)  ♦  DEL2 ( I ) 


SIMMER  LOOP 


Call  RAD  I  US  (  laMDA  ,  RANGE .  3) 


!*/G*  (AVC-PWR  i  TRAN  i  IREL-.NI);/(PI  i  R  i  R) 
iLAS£R=  AVGPWRAPI  *  R0  *  R0) 

?RGPEFF=  IAVG/ILASER 

IF  (  PROPEFF  .GT.  PBEST  )  THEN 
PBEST-PROPEFF 
vpbest=v 
END  IF 

CALL  FELOPT  (ETA.LAMDA) 

EFF=PROPEFF  i  ETA  *TGTEFF( LAUDA) 

IF  (EFF  .GT.  EBEST)  THEN 
EBEST=EFF 
VBEST=V 
END  IF 

WRITE  (19)  V.EFF 

IF  (IFLAG  .NE.  2)  THEN 
GO  TO  128 
END  IF 

*END  OF  WAVENUMBER  LOOP 

380  WRITE(6,390)  VPBEST.PBEST 
398  FORMAT ( '  BEST  WAVENUMBER  FOR  PROPAGATION', F10. 3, 
1  '  WITH  EFFICIENCY  OF  '.E12.4) 

WRITE  (6.400)  VBEST, EBEST 
400  FORMAT  ('  BEST  OVERALL  WAVENUMBER*' ,F 10. 3. 

1  '  WITH  EFFICIENCY  OF  '.E12.4) 


mimnmmmummmmmmmmmmmrimmm 

SUBROUTINE  MDLATM 

*  THIS  SUBROUTINE  READS  IN  THE  ATMOSPHERIC  PROFILE  FROM  'MODEL' 

REAL  P(33) ,T(33) ,ALT(33) ,P8Afi(32) .TBAR(32) ,RH0BAR(32) 

COrtlCN/ATM/P . T , ALT . PBAR . TBAfi , RHOBAR 

READ  (15)  (  P(I) ,T(I) ,ALT(I) ,  1=1.33) 

READ  (15)  (PBAR(I) ,TBAR(I) .RHOBAR(I) .  1=1.32) 

*  CONVERT  ALTITUDE  TO  METERS  AND  DENSITY  TO  KG/M3 


WSWiwWW 


DO  10  M=i,32 
riLT(N>=ALT(N)  i  1006.0 
RHG3AR ( N) =RHOBAR (N)  i  1000.0 
10  CONTINUE 

ALT ( 33) =HLT ( 33>  *  1080.0 

RETURN 

END 

SUBROUTINE  FCREAD  (V.  IFLAG, U1.U2) 

i  THIS  SUBROUTINE  READS  IN  MOLECULAR  ABSORPTION  DATA 
i  FROM  THE  FILE  GENERATED  BY  FASCODE.  IF  ENO  OF  FILE  IS 
*  REACHED.  IFLAG  IS  SET  TO  ONE. 

REAL  V,U1,V2,M0LABS(32) .AEXT(32) ,ABBS(32) ,RSCAT(32) ,0Z0NE(32) , 
1  C0NT(32) 

INTEGER  IFLAG, IENTER 

COMMON/OPDPTH/MOLABS  ,AEXT ,  ABBS ,  RSCAT ,  02  ONE .  CONT 
DATA  I ENTER/8/ 

19  READ  ( 14 ,END=  189)  V.(MOLABSU) .  1=1.32) 

IF  (IENTER  .EQ.  0)  THEN 
IF  (  0  .GE.  VI  )  THEN 
WRITE  (4.50)  V 

50  FORMAT ( '  STARTING  WAVENUMBER  WILL  BE  '.F10.4) 

ELSE 
GO  TO  10 
END  IF 
ENO  IF 

IENTER=1 


IF  (  V  .LE.  AINT(V2)  )  THEN 
RETURN 
ENO  IF 


108  WRITE  (4,150)  U 

150  FCRl-AT  C  LAST  l<W)ENUM8ER  READ  BY  FCREAD=  ',F16.- 
IFLAG®1 
RETURN 
END 

SUBROUTINE  LTREAD  (V, IFLAG) 

i  THIS  SUBROUTINE  REAOS  IN  ABSORPTION  WD  SCATTERING  DATA 
*  CREATED  BY  L0WTRAN4  AND  PERFORMS  INTERPOLATION 
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3LAL  V.Vl.VZ.RSCAT'.SZ)  .AtxT(  32.'  .a86S( 32)  .0ZGhE'.32>  .CENT 

1  ,RSCAi2(32)  .AtX72(32)  ,h8BS2(32)  ,0ZGNE2(32;  .CGNT2-' 32.! 

2  .RSCATK32)  .hE<Th32)  ,h8SS1(32)  .02CNEH32)  .CENTK32: 

3  .N0LABSC32) 

INTEGER  IEHTER, IFLrtG 

COtlON/OPOPTtVMOLABS.AEXT  .ABBS.RSCAT ,  OZONE.  CQNT 
DATA  IENTER/8/ 

IF  (I ENTER  .EQ.  8)  THEN 

READ ( 3) V 1 , < RSCAT  1(1)  .AEXTKI)  .ABBS  1(1)  .OZONE  1(1) . 

1  CQNT  1(1)  ,1=1,32) 

WRITER,  180)  VI 

188  FORMAT ( '  FIRST  WAVENUMBER  READ  BY  LTREAD=' ,F18.4> 
READ(3)V2,(RSCAT2(I) ,AEXT2(I) ,ABBS2(I) ,0ZGNE2(I) . 

1  C0NT2(I) ,1=1,32) 


I ENTER3 1 

IF(V  .LT.  VI)  THEN 

WRITE  <4,'(*  LTREAD  ERROR  *) ') 
STOP  2080 
ENO  IF 
END  IF 


150  IF  (V  .GT.  V2)  THEN 
V1=V2 

DO  280  1=1.32 
RSCAT1(I)=RSCAT2(I) 

AEXT1(I)=AEXT2(I) 

ABBS1(I)=ABBS2(I) 

OZONE  1(I)=OZONE2(I) 

C0NT1(I)=C0NT2(I) 

200  CONTINUE 

REAO<3.END=508)  V2.(RSCAT2(I)  .AEXT2U)  .ABBS2U)  ,OZONE2(I) . 
1  CONT(I) .1=1,32) 

GO  TO  150 

END  IF 

FACT4=<  VH4.0  -  V1H4.0  )/(  V23f*4.0  -  V1H4.8  ; 
FACT0R=(V-V1)/(V2-V1) 

DO  300  1=1,32 

RSCAT(I)=RSCAT1(I)  ♦  FACT4  {  (RSCAT2U)  -  RSCATKD) 
AEXT(I)=AEXT1(I)  +  FACTOR  *  (AEXT2U)  -  AEXTKI)) 

ABBS(I)  “ABBS  KI)  +  FACTOR  *  (A6BS2(I>  -  ABBS  1(D) 
OZONE(I)=OZONE1(I)  +  FACTOR  i  (0Z0NE2(I)  -  OZONEKD) 
CONTI I)=CCNT HI)  ♦  FACTOR  \  <X0NT2(I>  -  CONTKD) 


300  CONTINUE 


RETURN 


508  WRITER, 408)  U2 

680  FORMAT ( '  LAST  RENUMBER  READ  BY  LTREAD=- .F18.4) 

IFLAG=2 

RETURN 

END 

mmmm**m**m**mmm*mmmm*m*m*mm* 

SUBROUTINE  RADIUS  ( LAMDA, Z.R) 

i  THIS  SUBROUTINE  COMPUTES  THE  1/E2  RADIUS  OF  THE  BEtfl  AT 
l  A  SPECIFIED  DISTANCE  FROM  THE  LASER 

REAL  LAMDA . HOI , J ITTER , 3EAMQ , CGMPEN , FOCRNG , 0 1  ST , RHOTURB 

Cffl+ION/RAD/HOL  .JITTER,  BEAMQ .  FOCRNG ,  RHOTURB .  R0 .  COMPEN 


PI=3.141a 

ROGT2=SQRT(2.0) 

DlST=Z-HOL 

JITTER=JITTER  *  1.0E-4 

i  IF  FOCAL  RANGE  IS  ZERO,  THE  BEAM  IS  FOCUSED  AT  INFINITY 

IF  (ABS(FOCWG-0.0)  .LE.  0.91)  THEN 
FOCRNG=1.0E+25 
END  IF 

i  COMPUTE  TURBULENCE  RADIUS 

RTS= ( ( R00T2*DI ST*  ( 1 . 0-CQMPEN) *  LAMDA) /(PI* RHOTURB) )  3H2 . 8 

i  COMPUTE  JITTER  RADIUS 

RJS=  2.0*  (  (JITTERJfDIST)  **  2.0  ) 

*  COMPUTE  DIFFRACTION  AND  FOCUS  EFFECTS 

RDS=( (BEAMQ  *  LAMDA  *  DIST)/(PI  *  R0  *  R0)»  **  2.0 
RDS=  RDS  +  (  1.8-  (OIST/FOCRNG) )  **  2.0 
RDS=  RDS  *  R0  *  R0 

i  USE  ROOT  SUttt  SQUARE  TEQfUQUE  TO  COMPUTE  FINAL  RADIUS 
R=  SORT (  RDS  +  RJS  ♦  RTS  ) 

RETURN 

END 

*********************************************************** 
REAL  FUNCTION  REFIDX  (P,T,  LAMDA) 


‘  '-IS  FUNCTION  COMPUTES  THE  INDEX  OF  REFRACTION 


PEAi_  P.T.lAMOA 


RE^!DX=(5.dSE-l'?  <  ?V(T  i  lAMDa  s  lAMQA* 

RE? I DX=  REF I DX  ♦  (7.7cE-5  *  ?  /  T  •  +  1.8 

END 

REAL  FUNCTION  ONOT  (P.T.WHDA) 

*  THIS  FUNCTION  COMPUTES  THE  CWW6E  IN  REFRACTIVE  INDEX 
i  FOR  A  CHANGE  IN  TEMPERATURE 

REAL  P,T, LAMDA 

DNDT=-(5.35E-19  *  P)/(T  *  T  i  LAMDA  *  LAMDA) 

ONDT=ONDT  -  ((7.76E-5  i  P)/(T  i  T)) 

END 

REAL  FUNCTION  IREL  TNI) 

i  THIS  FUNCTION  COMPUTES  THE  FRACTION  OF  POWER  LEFT  IN  THE  BEAM 

*  AFTER  THERMl  BLOOMING  MS  TAKEN  AFFECT.  THE  INPUT  PARAMETER 
■*  IS  GESHARDT'S  DISTORTION  NUMBER 

REAL  NI 

IF  (NI  .LE.  8.1)  THEN 
I REL= 1.0 
RETURN 
END  IF 

IF  (NI  .GE.  40.8)  THEN 
IREL=0.8 
RETURN 
END  IF 

IREL3  1. 8/(1. 8  +  3.0625  X  NI  X  NI) 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
REAL  FUNCTION  VX(Z) 

X  DETERMINES  THE  HORIZONTAL  WIND  AT  ALTITUOE  Z 
REAL  Z 


UX=8.0 


-  «  -*»  —  ,  1  ■**  -  •  ‘ 


.^tLvvV  lV-..--vO vC" s-  vlvl  .C'  +. 
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VX=27.2  -  2. 72E-4  i  Z 
RETURN 

ELSE  IF  (Z  .GE.  12288.3)  THEN 
UX=82.4  -  2.33E-3  *  2 
RETURN 
ELSE 

VX=8.75  +  3.42E-3  I  Z 
RETURN 

END  IF 
END 

SUBROUTINE  FELOPT(ETA,UWOA) 

6  THIS  SUBROUTINE  COMPUTES  THE  DEVICE  EFFICIENCY  FOR 
*  A  GIVEN  WAVELENGTH 

SEAL  ETA.LAMOA.LAWAP.NjIPEAKB.IBARM.IB.EPSLON.THETAM, 

1  IN, PBAfto.PBARL.GAWIAB.GAlflAR, KAPPA. X.G.FK, DELTA, VALUE 

if  DEFINE  CONSTANTS  AND  NOMINAL  DEVICE  VALUES 

PBARA=2.1E+4 

PBARL=3.8E+4 

LAM0AP=9.8E-3 

N=2.8 

IBARM=1.8E+11 
IPEAKEN  158.0 
10=1354.0 
LM=50.8 
THETAM=1.0£-4 
EPSLd+=0 . 1 
PI=3.1414 
GAM1AB=10.0 

ifCQMPUTE  KAPPA 

KAPPA=(  ( 4 . 03(P!  ifEPSLON)  /  ( LMJfLMiLAMDAifTHETAM) )  U .  447 
KAPPA=KAPPA3f(LAM0AP/t4.03fPI3fN))  36X2.0 

if  DETERMINE  FROM  CURVE  FITS  OPTIMUM  VALUES  OF  X.G.FK  GIVEN  KAPPA 

X=( -8 .852836 KAPPA* KAPPA)  ♦  <  4 . 828436 KAPPA)  ♦  .52725 

6=  <  -3 . 449536 KAPPASf KAPPA)  ♦  ( 1.4418i(KAPPA)  ♦  1.74241 

FK=((EXP(-G/X)  36  X  36  X)  +  KAPPA)  3636  (-1.5) 

FK=FK  6  (EXP(-G/X)  6  X  6  X) 


FK-FK  *  (  uG*G.>  -  1.0)/(G*G)  >  *  \  1  !•  a.)  xi  1.5  1 
FK=FK  i  (  KAPPA  **  1.25  ) 


*  COMPUTE  HIGGLER  VECTOR  POTENTIAL 

PW=(KAPPA**(-.5))  *  X  k  EXP(-G/(2.0StX)> 

J COMPUTE  RESONANCE  GtftIA  AT  ENTRANCE  TO  WIGGLER 

GM»R=(  (1.0+ <  AW*#1) )  *  LA1DAP*X)  /  ( UW)A*2 . 0*  ( KAPPA**  ( .  5) ) ) 
GAMMAR=SQRT  (GAWAR) 

iCOMPUTE  DELTA 

DELTA=(  (LAMOAP/LAMDA)  **  1 .5)  *<GW+ttB- 1 .0)  *PBARL*I0 
DELTA=DELTA/  ( LN*LM*  I  BARM*  I  PEAKBiFK*  ( GAftttR- 1.0)) 

*  FIND  VALUE  FROM  CURVE  FIT 

VALUE=(ll.9642*DELTA*DELTA)  -  (3.7919*0ELTA)  +.55234 

*  COMPUTE  EFFICIENCY 

ETA=((  1.0/UALUB  +  (PBARA/PfiARU)**(-1.0) 

RETURN 

END 

************************************************************* 
REAL  FUNCTION  TGTEFF(LAMDA) 

*  THIS  FUNCTION  SCALES  RELATIVE  EFFICIENCY  OF  TARGET  ABSORPTION 

*  WITH  .3448  MICRONS  BEING  THE  BEST 

REAL  LAMDA 


TGTEFF=SQRT ( .3448E-6/LAMDA) 
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